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I.       INTRODUCTION 

A.  GENERAL 

The  purpose  of  this  study  was  to  determine  the  feasibility  of 
transporting  the  Ml  19  howitzer  in  the  firing  position  over  terrain 
and  at  speeds  used  during  normal  tactical  operations.  The  end 
objective  was  to  determine  if  costly  design  changes  would  be 
required  to  safely  permit  towing  the  Ml  19  Howitzer  in  the  firing 
position.  The  approach  used  was  to  model  the  vertical  plane 
dynamics  of  the  towing  vehicle  (HUMMV)  and  the  howitzer,  model 
representative  terrain  inputs  to  the  wheels  in  the  form  of  potholes, 
and  to  perform  experiments  to  validate  the  model.  It  is  proven  that 
the  simulation  model  is  reasonable  and  it  follows  that  stresses  in  the 
critical    components    are    not    expected    to    warrant    redesign. 

B.  BACKGROUND 

The  Ml  19  Howitzer  is  a  British  designed  howitzer  that  is 
currently  being  built  at  Rock  Island  Arsenal  (RIA),  Illinois  and 
fielded  to  US  Army  Field  Artillery  units.  The  Ml  19  Howitzer 
operational  manuals  specify  towing  procedures  that  require  the  gun 
tube  to  be  placed  in  the  stowed  position  when  moving  the  howitzer 
during  tactical  operations.  These  towing  procedures  are  different 
and  more  difficult  to  implement  than  those  used  for  all  other  towed 
howitzers  in  the  US  Army.  All  other  US  Army  towed  howitzers  are 
towed     during     tactical     operations     with     the     gun     tube     in     the     firing 
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position,  facilitating  rapid  emplacements  and  displacements  that  are 
required    for    successful    tactical    operations. 

The  Ml  19  Howitzer  has  an  existing  A-Frame  brace  which  can 
support  the  gun  tube  when  the  howitzer  is  towed  in  the  firing 
position.  The  primary  concern  of  this  study  was  whether  or  not  the 
A-Frame  and  the  T-Bar  structure  on  the  howitzer's  trail  section,  to 
which  the  bottom  of  the  A-Frame  brace  connects,  will  be  strong 
enough  to  support  the  loads  generated  when  towing  the  howitzer 
over  rough  terrain.  This  study  also  addresses  concerns  about  trail 
loads  on  the  howitzer.  This  study  does  not  evaluate  rollover  stability 
for   the   Ml  19   Howitzer   when    towing    in    the   firing   position. 

The  terrain  that  the  HUMMV/Howitzer  system  can  traverse  is 
clearly  limited  by  the  terrain  the  HUMMV  can  negotiate,  therefore, 
the  Ml  19  Howitzer  and  HUMMV  have  been  evaluated  as  a  system.  A 
dynamic  simulation  model  of  this  system  was  seen  to  be  a  necessary 
tool  to  aid  in  making  design  decisions  from  a  rational  viewpoint  and 
is   the  basis   of  this   study. 

C.     METHOD  OF  ANALYSIS 

To  analyze  the  Ml  19  Howitzer  a  computer  model  of  the 
HUMMV/Howitzer  system  was  developed  to  predict  the  response  of 
the  system  to  roadway  inputs.  This  computer  model  treats  the  major 
mass  elements  of  the  HUMMV/Howitzer  system  as  lumped  masses  to 
develop  vertical  plane  heave  and  pitch  equations  of  motion  for  the 
system.  The  development  of  this  computer  model  is  contained  in 
Chapter    II.       To    drive    the    computer    simulation    model,    terrain    models 


were  developed  for  both  a  "worst  case"  trench  type  pothole  and  a 
sinusoidal  "washboard"  roadway.  Details  on  the  development  of 
these  roadway  inputs  are  contained  in  Chapter  III.  Chapter  IV 
contains  the  results  of  driving  the  computer  model  over  the  roadway 
inputs.  Chapter  V  explains  the  structural  analysis  that  was  used  to 
determine  the  maximum  stress  in  the  A-Frame  and  the  T-Bar. 
Chapter  V  also  contains  the  procedures  and  results  of  calibrating  an 
A-Frame  that  was  used  as  a  load  cell  for  an  instrumented  dynamic 
test.  Chapter  VI  explains  the  conduct  and  the  results  of  the 
instrumented  test  on  the  Ml  19  Howitzer/HUMMV  system  that  was 
conducted  at  Aberdeen  Proving  Grounds,  Maryland  in  April  1993. 
Also,  in  Chapter  VI  the  results  of  the  model  are  compared  to  the 
experimental  results  from  the  instrumented  test.  This  comparison 
showed  reasonable  correlation  between  the  model  and  experimental 
results.  Conclusions  and  recommendations  based  on  this  study  are 
contained  Chapter  VII.  The  final  recommendation  is  that  towing  in 
the    firing    position    should    be    allowed. 


II.       IIUMMV/HOWITZER    SYSTEM    COMPUTER   MODEL 

The  computer  model  developed  is  based  only  on  the  heave  and 
pitch  equations  of  motion  (EOM)  for  the  HUMMV/Howitzer  system. 
The  reason  for  this  simplification  is  that  the  lateral,  yaw,  and  roll 
motions  are  much  less  significant  for  determining  the  forces  in  the 
A-Frame  for  the  type  terrain  inputs  considered  by  this  study.  While 
horizontal  plane  motions  are  key  elements  for  the  analysis  of  turns 
and  rolling  during  transit,  the  concerns  of  this  study  indicate  that  a 
vertical  plane  model  is  adequate.  The  major  mass  groups  of  the 
HUMMV/Howitzer  were  treated  as  lumped  masses  in  the 
development  of  the  equations  of  motion.  The  equations  of  motion 
were    developed    using    Lagrange    techniques. 

A.       HUMMV/HOWITZER    SYSTEM    PARAMETERS 

The  primary  sources  of  information  for  system  parameters  for 
the  HUMMV  and  howitzer  were  the  US  Army  Tank-Automotive 
Command  Research,  Development  &  Engineering  Center  Technical 
Report  No.  13337,  [Ref.  1],  for  HUMMV  data  and  US  Army  Combat 
Systems  Test  Activity,  APG,  Report  No.  P-83179,  [Ref.  2],  for  the 
Ml  19  Howitzer  data.  Additionally,  many  system  parameters  were 
determined  from  direct  measurement  on  the  HUMMV  and  howitzer 
or  from  design  diagrams,  or  obtained  from  the  study  sponsor  at  Rock 
Island  Arsenal  (RIA).  For  readers  of  this  study  that  are  unfamiliar 
with    the    components    of    the    Ml  19    Howitzer,    Appendix     G     contains 


several  diagrams  that  depict  the  howitzer  components  that  are 
important    for    this    study. 

The  system  parameters  used  for  the  HUMMV  are  based  on  using 
the  Model  No.  M1069  HUMMV  as  the  prime  mover  for  the  Ml  19 
Howitzer.  The  cargo  loading  of  the  HUMMV  significantly  effects  its 
center  of  gravity  and  consequently  effects  all  of  the  HUMMV's 
system  parameters  which  are  dependent  on  its  center  of  gravity. 
Modeling  was  therefore  limited  to  using  the  same  cargo  loading 
configuration  that  was  used  in  [Ref.  1]  and  other  loading  effects  were 
ignored.  These  HUMMV  simplifications  are  justified  because  for  this 
computer  model  the  HUMMV  serves  primarily  as  the  driving  force  at 
the  howitzer  hitch  causing  pitching  motions  of  the  entire  howitzer 
and  we  are  not  concerned  with  the  loading  of  HUMMV  components. 
The  major  components  of  the  HUMMV/Howitzer  system  computer 
model    are    identified    in    Figure    2-1    below. 

In  the  model,  the  tires  act  as  linear  springs  which  provide  force 
inputs  to  the  HUMMV  and  howitzer  wheels,  based  on  the  terrain  over 
which  the  model  was  "driven".  The  suspension  systems  for  both  the 
HUMMV  and  the  howitzer  are  modeled  as  springs  and  dampers 
connected  in  parallel.  The  A-Frame  and  T-Bar  are  modeled  as  a  set 
of    stiff    springs    connected      in      series.  The      following     paragraphs 

provide  more  detail  on  the  modeling  of  the  individual  components 
and  the  determination  of  appropriate  parameter  values.  Appendix  A 
provides  a  listing  of  the  numerical  values,  units,  and  the  symbols 
used     for     each     system     parameter     in     the     computer     model.  The 

parameter    values    for    those    components    with    identical    components    on 


Figure    2-1    HUMMV/Howitzer    Model    Components 


the  right  and  left  sides  of  the  system  (wheels  and  suspension 
components)  are  multiplied  by  two.  The  values  in  Appendix  A  for 
these    components    reflect    the    effects    of   both    sides. 

1.  HUMMV    and    Howitzer    Tires 

The  HUMMV  has  the  same  tires  on  the  front  and  rear, 
however,  the  front  and  rear  tire  pressure  is  different  (20  psi  vs  30 
psi)  resulting  in  the  rear  tire  spring  being  much  stiffer.  The 
stiffness  used  for  the  HUMMV  tires  was  obtained  from  normal  force 
vs  tire  deflection  data  [Ref.  1].  The  stiffness  of  the  howitzer  tire  used 
in  the  model  was  estimated  because  no  measured  data  on  the 
howitzer    tire    stiffness    was    available. 

The  tires  are  modeled  as  linear  springs  which  transmit  a 
vertical  force  to  the  respective  wheel  masses  proportional  to  the 
amount  that  the  springs  are  stretched  or  compressed.  The  amount  of 
force  is  determined  by  the  relative  dynamic  displacement  between 
the  terrain  input  for  the  wheel  and  the  center  of  mass  of  the  wheel. 
If  the  tire  leaves  contact  with  the  ground  there  will  be  no  input  from 
the  terrain  until  it  again  returns  to  contact.  To  account  for  this  loss 
of  contact  situation  in  the  computer  model,  the  tire  spring  force  is  set 
to  the  static  weight  that  the  tire  supports  when  the  tensile  (de- 
compressive) spring  force  exceeds  the  static  weight  on  the  tire.  The 
tire  springs  are  turned  back  on  in  the  computer  model  when  the 
conditions    are    met    for    road    surface    contact    to    be    re-established. 

2.  HUMMV     Suspension     System 

The  HUMMV  suspension  system  consists  of  a  shock 
absorber    and    coil    spring    for    each    wheel.       The    shock    absorber   passes 


through  the  middle  of  the  coil  spring.  The  force  in  the  spring  is 
linearly  proportional  to  the  relative  displacement  between  the  center 
of  mass  of  the  wheel  and  the  attachment  point  on  the  HUMMV  body. 
The  values  for  the  spring  constants  were  obtained  from  Table  5-13 
of  [Ref.  1].  Precise  modeling  of  the  shock  absorbers  is  difficult 
because  they  have  a  mid-stroke  region  and  a  hydraulic  bump  stop 
region.  In  the  mid-stroke  region  the  shock  absorber  force  is  linearly 
dependent  on  the  relative  velocity  of  the  center  of  mass  of  the  wheel 
and  the  attachment  point  on  the  HUMMV.  However,  the  linear 
damping  constant  changes  depending  on  whether  the  force  in  the 
shock  absorber  is  tensile  or  compressive.  When  the  shock  absorber 
is  compressed  or  stretched  into  the  hydraulic  bump  stop  region  the 
damping  of  the  shock  absorber  greatly  increases.  This  sudden 
increase  in  damping  causes  problems  in  the  integration  routine  that 
was  used  for  the  computer  model.  To  overcome  these  problems  and 
still  account  for  the  energy  loss  due  to  the  increased  damping  in  the 
hydraulic  bump  stop  region,  the  mid-stroke  linear  damping  constant 
obtained  from  [Ref.  1]  was  increased  by  20%  for  both  the  front  and 
rear  shock  HUMMV  absorbers.  Figures  2-2  and  2-3  shows  the  force 
versus  velocity  plots  that  were  obtained  from  [Ref.  1]  for  the  front 
and  rear  shock  absorbers,  respectively. 
3.         Howitzer     Suspension     System 

Modeling  the  howitzer  suspension  system  is  more 
complicated  than  the  HUMMV  because  it  consists  of  a  torsional  spring 
and  a  shock  absorber  that  is  not  vertically  oriented.  The  torsional 
spring    is    modeled   as   a   linear   vertical    spring    that   is   connected   to   the 
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Figure  2-2     Force  vs   Velocity   Plot   for  the  Front  HUMMV   Shock 
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Figure  2-3      Force  vs  Velocity  Plot  for  the  Rear  HUMMV   Shock 
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center  of  mass  of  the  howitzer  wheel  and  attached  to  the  howitzer 
trail  directly  above  the  wheel  center  of  mass.  The  effective  spring 
stiffness  constant  was  calculated  based  on  the  measured  data  of  the 
torsional  spring  stiffness.  Figures  2-4  through  2-7  show  the  angular 
deflection  versus  vertical  force  plots  that  were  contained  in  [Ref.  2]. 
Using  the  average  torsional  spring  stiffness  from  this  data  and  the 
length  of  the  road  arm,  an  angular  rotation  was  converted  to  an 
equivalent  vertical  displacement  at  the  end  of  the  road  arm.  The 
spring  constant  for  the  model  was  then  calculated  by  determining 
the  constant  value  that  must  be  multiplied  times  this  vertical 
displacement  at  the  end  of  the  road  arm,  to  produce  a  moment  that  is 
equivalent  to  that  caused  by  the  torsional  spring  for  a  specified 
angular     rotation. 

A  similar  procedure  was  used  to  model  the  howitzer's 
shock  absorber  as  a  vertical  shock  absorber  with  a  linear  damping 
constant  that  is  velocity  dependent.  The  attachment  points  for  the 
ends  of  the  shock  absorbers  in  the  model  are  the  same  as  the 
modeled  end  points  of  the  howitzer  suspension  spring.  The  actual 
non-vertical  shock  absorber  force  was  transformed  to  a  vertical  force 
in  the  model  by  determining  the  amount  of  vertical  force  applied  a  t 
the  end  of  the  road  arm  needed  to  produce  an  equivalent  moment 
about  the  road  arm  pivot  point.  Like  the  HUMMV  shock  absorbers, 
the  howitzer  shock  absorbers  are  non-linear  and  have  different 
damping  values  for  tension  and  compression  that  are  accounted  for 
in    the    computer     model.  This     non-linear     damping     in     the     shock 

absorbers    is    a    significant    aspect    of   the    computer    model.       The     values 
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used  for  the  shock  absorber  damping  constants  were  determined 
from  force  versus  velocity  test  data  provided  by  RIA.  Figure  2-8 
shows    the    force    versus    velocity    plot    for    the    howitzer    shock    absorber. 
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Figure     2-4     Torsional     Spring     Stiffness     for     Serial     No.     22     (Upward 
Motion) 
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Figure    2-5       Torsional    Spring    Stiffness    for    Serial    No.     22     (Downward 
Motion) 
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Figure     2-6    Torsional     Spring     Stiffness     for     Serial     No.     27     (Upward 
Motion) 
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Figure    2-7    Torsional     Spring     Stiffness    for    Serial    No.     27    (Downward 
Motion) 
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4.  A-Frame/T-Bar      Stiffness 

The  A-Frame  and  T-Bar  are  modeled  as  stiff  springs 
connected  in  series  to  determine  a  combined  linear  spring  stiffness. 
The  stiffness  of  the  A-Frame  was  determined  during  the  A-Frame 
calibration  test  (explained  in  Chapter  V).  The  stiffness  of  the  T-Bar 
was  determined  by  using  the  principle  of  superposition  of  forces  to 
determine  the  effective  bending  stiffness  for  the  non-symmetric 
bending  of  the  T-Bar  that  is  caused  by  the  A-Frame.  This  non- 
symmetric  bending  of  the  T-Bar  results  because  the  angle  at  which 
the  A-Frame  transmits  forces  and  moments  to  the  T-Bar  is  not 
perpendicular    to    an    axis    of   symmetry    of   the    T-Bar. 

5.  Equilibrator     Springs 

The  equilibrator  springs  provide  only  small  dynamic  forces 
because  their  combined  stiffness  is  nearly  1000  times  less  than  the 
effective      stiffness      of     the      A-Frame/T-Bar     assembly.  Therefore, 

equilibrator    springs    are    neglected    in    the    computer    model. 

6.  Elevating    Clutch 

The  effects  of  the  elevating  clutch  of  the  howitzer  were 
neglected  in  the  final  computer  model.  It  is  believed  backlash  in  the 
elevating  mechanism  gears  will  prevent  the  clutch  from  being 
activated  while  towing  with  the  A-Frame  attached  because  the 
relative  rotation  of  the  gun  tube  is  very  small.  However,  because  of 
some  uncertainty  about  the  clutch  several  computer  runs  were  made 
with  clutch  modeled  as  a  linear  torsional  damper.  The  value  of  the 
torsional  damper  used  in  the  model  was  estimated  because  no  test 
data    on    the    clutch    was    available.       These    model   runs    with    the    clutch 
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included,  produced  only  slightly  different  results  from  the  runs 
without  the  clutch  included.  Therefore,  because  of  the  previously 
mentioned  belief  that  the  clutch  would  not  be  activated,  it  was 
neglected    in    the    final    model. 

7.  HUMMV     Component     Mass    Properties 

The  HUMMV  is  divided  into  three  lumped  masses,  the 
HUMMV  body  mass  and  two  wheel  masses.  The  front  and  back 
HUMMV  wheels  are  treated  as  having  identical  mass  values  in  the 
model.  Each  wheel  weight/mass  includes  the  weight  of  the  tire, 
wheel,  shocks,  suspension  springs,  and  suspension  support  arms.  The 
HUMMV  body  weight/mass  varies  dependent  upon  it's  cargo  load. 
The  location  of  the  center  of  gravity  also  can  vary  and  is  determined 
by  the  positioning  of  the  cargo.  Therefore,  the  values  used  in  the 
computer  model  for  the  HUMMV  mass,  the  mass  moment  of  inertia 
for  pitch,  and  lengths  that  are  based  dependent  on  the  center  of 
gravity,  are  estimated  values.  These  estimated  values  are  based  on 
the  same  loading  situation  as  was  used  in  [Ref.  1].  The  pitch  mass 
moment  of  inertia  for  the  HUMMV  is  about  the  center  of  gravity  of 
the  sprung  HUMMV   mass. 

8.  Howitzer     Component     Mass     Properties 

The  howitzer  is  divided  into  three  separate  lumped  masses 
in  the  model;  the  howitzer  wheel  masses,  the  elevating  mass,  and 
the  combined  trail/saddle  mass.  The  wheel  weight/mass  includes 
the  tire,  rim,  and  approximately  30%  of  the  road  arm  that  rests  on 
the  wheel  assembly.  The  elevating  mass  includes  the  gun  tube,  the 
cradle,   the  recoil  system    and    all    other    parts    that    move    the    when    the 
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gun  is  elevated  about  the  trunion.  The  mass  moment  of  inertia  of  the 
gun  tube  for  rotation  about  the  trunion  was  estimated  using  the 
approximate  masses  of  all  the  elevating  components.  The  center  of 
gravity  location  for  the  elevating  mass  was  estimated  by  RIA.  The 
combined  trail/saddle  mass  includes  all  of  the  gun  mass  except  the 
elevating  mass  and  the  wheel  mass.  The  center  of  gravity  location 
for  the  trail/saddle  mass  and  the  pitching  mass  moment  of  inertia  for 
the  trail/saddle  mass  was  estimated  from  design  diagrams. 
9.         Length     Properties 

All  lengths  used  in  the  computer  model  are  identified  in 
Figure  2-9  below.  These  length  properties  were  determined  either 
by  direct  measurements  on  the  howitzer  or  by  scaling  from  design 
diagrams. 

B.      EQUATIONS   OF  MOTION   (EOM) 

In  this  section  the  heave  and  pitch  equations  of  motion  (EOM) 
for  the  model  are  developed  from  Lagrange's  equations  using  the 
HUMMV/Howitzer  system  parameters  that  were  previously 
explained.  Lagrange's  equations  require  the  determination  of  the 
system's  total  kinetic  and  potential  energy  and  the  power  dissipated 
by   the   dampers,    in    terms   of   a    set   of   generalized   coordinates. 

Throughout  this  development  of  the  EOM,  expressions  are 
developed  that  are  functions  of  generalized  positions  and  velocities, 
angular  displacements,  angular  velocities,  and  angular  accelerations. 
To  simplify  the  EOM  with  minimal  loss  of  representational  accuracy, 
these    expressions    are    linearized    using    small    angle    approximations. 
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Figure    2-9    Length    Properties    of    HUMMV/Howitzer    System 
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The  sine  of  a  small  angle  is  approximated  as  just  the  angle  (sin0=0) 
and  the  cosine  of  a  small  angle  is  approximated  as  equal  to  one 
(cos6=l).  These  approximations  are  valid  because  the  pitch  or 
angular  rotations  are  predicted  to  be  less  than  five  degrees.  Order  of 
magnitude  considerations  are  used  to  further  simplify  the  EOM. 
Because  the  angular  displacements  are  small  their  angular  velocities 
and  angular  accelerations  are  also  small.  Therefore,  any  terms  that 
contain  an  angular  displacement,  angular  velocity,  or  angular 
acceleration  multiplied  times  another  angular  displacement,  angular 
velocity,  or  angular  acceleration  are  neglected  because  these  terms 
will    be   very    small. 

The  general  form  of  Lagrange's  equations  used  is: 


d^ 
dt 


dT       au        dF_ 
dq{       dq{       dqi 


dT 

aq, 


where: 

qi    =    generalized    coordinates 

T   =   total   kinetic   energy   of  the   system 

U  =  change  in  the  potential  energy  of  the  system    with    respect    to 

its  potential   energy   in   the   static-equilibrium   position 
F   =   the   power   dissipated   by    forces   in    dampers 
Q   =  the     input  force  acting  on   the  system. 

The  development  of  the  above  components  of  Lagrange's  equation 
and  the  assembly  of  these  components  into  the  system's  heave  and 
pitch  EOM  follows.  Due  to  the  length  of  this  development  the  reader 
is  provided  only  the  final  results  from  each  step.  The  intermediate 
multiplication's,    differentiation's,    and    linearization's    are    omitted. 


2  1 


1.         Generalized     Coordinates     (qj) 

The  set  of  generalized  coordinates,  q;9  that  was  used  to 
uniquely  specify,  the  motion  of  the  HUMMV/Howitzer  system  are 
depicted    on    Figure    2-10    and    are    defined    follows: 

*  Z1&Z2*     The    vertical      displacements      for      the      front      and      rear 

HUMMV    wheel     masses,    respectively. 

*  Z3:      The  vertical   displacement     of  the  center  of  gravity   of  the 

HUMMV    body. 

*  Z4:      The  vertical   displacements     of  the  center  of  gravity  of  the 

howitzer    wheel    mass. 

*  9i :   The  pitch  angle  of  the  HUMMV  body  about  its  center  of 

gravity. 

*  82:      The   pitch   angle   of  the   howitzer   trail/saddle   mass   about   its 

center    of    gravity. 

*  63:      The  absolute  rotation   of  the  gun   tube. 

As  depicted  on  Figure  2-10,  upward  vertical  motions  and  clockwise 
rotational  motions  are  defined  as  positive.  The  equilibrium  position 
for  each  of  the  generalized  coordinates  is  the  static  equilibrium 
position  for  the  respective  components  that  are  correlated  to  each  of 
the  generalized  coordinates.  In  defining  positions  to  key  locations  in 
the  system  unit  vector  notation  of  i,  j,  and  k  will  be  used 
corresponding  to  the  x,  y,  and  z  axis  directions  of  the  Cartesian 
coordinate    reference    system    shown    on    Figure    2-10. 

The  kinetic  and  potential  energy  of  the  system  and  the 
damping    and    input    forces    are    all    computed    in    terms    of   these 
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generalized  coordinates  and  the  road  inputs  Ui,  U2,  and  U3.  To  carry 
out  the  required  computations,  position  vectors  to  key  locations  on 
the  HUMMV  and  howitzer  are  defined  in  terms  of  the  generalized 
coordinates.  The  position  vectors  to  these  key  locations  are  depicted 
on   Figure   2-11    and   defined   as    follows: 


* 


* 


n  &  r2    =    the    front   and   rear   attachment   points,    respectively,    for 
the    suspension    springs    and    shock    absorbers    on    the    HUMMV. 

r3=    the   center   of   mass   of  the   trail/saddle    mass 

r4=    the    modeled   attachment   location    of  the    suspension    spring 
and    shock    absorber    of   the    howitzer 


*  rs=   the   location   where   the   A-Frame   attaches   to   the   T-Bar 

*  r6=   the  location   of  the  center  of  mass  of  the  gun   tube/elevating 

mass 

*  rj=   the   location   where  the   A-Frame   attaches   to  the  gun 
tube/elevating      mass 

These  key  locations  that  are  defined  above  can  be  expressed  as 
position  vectors  in  terms  of  the  above  defined  lengths  (Figure  2-9), 
generalized  coordinates  (Figure  2-10),  and  the  unit  vectors  i,  j,  and  k 
of   the    Cartesian    coordinate    reference    system    (Figure    2-10)    as    follows: 

*  Tj  =  (Z3  -  L^inCSjJk  +  (l^cos^i 

*  r2  =(Zj  +  L2sin(61))k  -  (L2cos(0,))i 

*  r3  =(Z,  +  LjSintOj)  +   Lusin(62))k  -  (i^cos^)  +  Lncos(62))i 

*  r4  =(Zj  +  L3 sin^)  +   L12sin(02))k  -  (L3cos(01)  +  L12cos(02))i 

*  r5  =(Z3  +  LjSinfej  +  (L13+ L16)sin(02))k  -  (L, 008(0,)  +  (L,s  +  L, Jcos^))! 
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*  r6  =  (Z3  +  I^sin(0I)+L16sin(02)  +  LCGsin(0o  +  03))k  - 

*  r7  =  (Z3  +  l^sin(01)+L16sin(02)  +  LAsin(0o+03))k  - 


L3  008(0,) +  L,6cos(02)Y 

,+LCGcos(0o  +03)         / 
L3  cos(0[ )  +  L16cos(02)\ 
,  +  LAcos(0o+03)  / 


The  velocities  of  these  key  locations  can  be  determined  by 
taking  the  time  derivatives  of  the  position  vectors.  However,  only 
the  velocities  of  locations  where  there  are  velocity  dependent  actions 
occurring     in     Lagrange's     equations     need     to     be     determined.  The 

required    velocities    of   the    key    locations    on    the    HUMMV    and    howitzer 
are   as    follows: 


*  r,  =(z3  -  LjeosCOJoJk  -  (Lt sinCeje'Ji 

*  r2  =(z,  +  LjCOsfO^oJk  +  (L^in^)©,)! 

*r3=(Z3  +  Ljcos^)©!  +  Lncos(02)02)k  +  (l3  sin(Q1  )Qt  +  Lnsin(02)02)i 

*r4=(Z3  +  l^cosC©,)©!  +  L12cos(02)02)k  +  (i^sinC©,)©,  +  L12sin(02)02)i 


r,  = 


rZ3  +  L3cos(01)01  +  L16cos(02)02> 

+  ^008(00+03)0,  , 


k  + 


L3sin(01)01  +  L16sin(02)0 
+  LCGsin(0o  +  03)03 


i  \ 


In   addition   to   the   above   velocities,   the  vertical   velocity   of  the   CG   of 
the  HUMMV,   Z3,   and  each  of  the  wheel   masses,   Zp  Z2,  &  Z4    are 

needed    in    Lagrange's    equation. 
2.       Kinetic    Energy    (T) 

The    velocities    defined    above    are    used    to    calculate    the 
kinetic   energy   of  the    system.      The   kinetic   energy   is   calculated   by 
summing    the    linear   and    rotational    kinetic    energy    for    each    of   the 
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system's    masses.      The    total    kinetic   energy    of   the    system    is   given    by: 
I  T,=  £  (^m.v,2  +  fcl^2)         i  =  1,6;    j  =  1,3 

where: 

Tj  =  the    kinetic   energy    of  each    mass    in    the    system 

mi  =   the      mass   of  each   lumped   mass 

vj  =  the   linear   velocity   of  the   each    mass 

Ij  =  the    heave    mass    moment    of  inertia   of  each    rotating    mass 
about   its      center   of  gravity 

6j  =  the   angular   velocity    of   each    rotating    mass 

Therefore,  the  total  kinetic  energy  of  the  HUMMV/howitzer  system 
was  determined  by  summing  the  kinetic  energy  components  as 
follows: 

T  =  ^M1(Z1)2  +   ^M2(Z2)2  +  iM3(Z3)2+   ^M4(Z4)2+  ^M5(r3)2 
^M6(r6)2  +  il^e,)2  +  il2(02)2  +  il3(63)2 

In  the  above  expression  the  square  of  the  velocities  is  determined  by 
substituting  in  the  appropriate  velocity  expression  as  defined  above, 
and  taking  the  dot  product  of  the  velocity  with  itself.  This 
substitution  and  dot  product  operation  is  left  for  the  reader.  When 
taking  this  dot  product  all  of  the  k  dot  i  and  i  dot  k  cross  product 
terms  drop  out  because  the  i  and  k  unit  vectors  are  orthogonal. 
Having    determined    the    total    kinetic    energy,    the     time     derivative     of 

the    partial    derivative    of    the    total    kinetic    energy    with    respect    to    the 

d  dT 
velocity     of     each     of      the      generalized      coordinates,      -—     can     be 

dt  dqi 

determined.  Carrying  out  this  differentiation  and  linearizing  yields 
the    following: 
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d  dT 

dtaz,        '  ' 

d  dT 

* r-  =     M2Z2 

dtaz2       2  2 

--^-=  (M3+M5  +  M6)Z3  +  (M5L3  +  M6L,)0;  +  (M5Ln  +  M6L16)02 
*  dt  oZ3 


+   M6LCGcos(0o)03 


d   dT 

*  — ^-=   M4Z4 

dtaz4        4  4 

d  dT 


*  dt  ae 


r=   (M5L3  +  M6L3)Z3  +  (I,  +  M51V+    M6V)0, 


+  (M5L3Ln  +   M6L3L16)02  +   MjLjLcoCOsOoJBj 
tIt  =   ^Lii  +  M6L16)Z,  +  (M5L3Ln  +  M^Lje; 

*  dt  a62 

+  (I2   +  M5LU2+  M6L162)02  +  M6L16LCGcos(0o)03 

d  dT 

~ZJd~  =    M6LcGcos(0o)Z3  +  MgLaLcGCOsCOoje,  +    M6L16LCGcos(0o)02 

*  dt  a©3 

+  (I3  +  M6LCG2)03 

The  total  kinetic  energy,  T,  must  also  be  differentiated 
with  respect  to  each  of  the  generalized  positional  coordinates, 
however,  for  this  system  this  differentiation  is  equal  to  zero  for  all 
coordinates. 

3.        Potential     Energy     (U) 

The  potential  energy  for  this  system  results  from  the 
strain  energy  that  is  developed  in  the  suspension  springs  and  the  A- 
Frame  spring.  The  gravitational  potential  energy  that  arises  from 
changes  in  vertical  positions  is  ignored  because  the  resulting  terms 
drop  out  of  the  dynamic  equations  of  motion.  The  potential  energy 
in    any    spring    is    determined   as    follows: 
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U.^K.CAx,)2 

where: 

Ui  =    the    potential    strain    energy    in    each    spring 

Kj  =    the    spring    stiffness    for    spring 

Ax;    =   the  change   in   length   of  the   spring. 


Since  locations  to  the  attachment  points  of  the  ends  of  all  the  springs 
have  been  defined  above  in  a  common  inertial  coordinate  reference 
system,  taking  the  difference  of  the  position  vectors  of  the  end  points 
of  each  spring  will  yield  the  change  in  length  of  the  spring.  Thus  for 
this  system,  the  total  potential  energy  is  found  by  summing  the 
potential    energy    in   each    of  the   springs   as    follows: 


u4Ki(zrui)2  +  iK2(ri  -  zif  +  Kfe-^f  +  iK4(r2  -  zi) 

+  iK6(Z4-U3f +  7K7(r4  "  Z4f+  7KA(r7  -   *5f 


Therefore,  the    partial    derivative    of    the    total     potential     energy     with 

respect    to  each    of    the    generalized    coordinates,     — ,     after     linearizing 

results    in  the    following: 

eu 


* 


dZ, 


=    KjZj    +    K2Zj  -  K2Z3  +   K^LjSj  ~~  KjUj 


*  £L  =  K3Z2    +   K4Z2  -    K4Z3  -  K.L^  -   K3U2 

az2 

-— -    -K2Z,  -  K4Z2  +  (K2  +  K4  +  K7)Z3  -  K7Z4 

*  dZ3 

+  (K2Lj  +   K4L2  +  K7L3)0,  +   K  7L,202 

*£-Tm   "K7Z3  +  (K6  +  K7)Z4  -   K7L301  -  K7L1202  -   K6U3 

oZ4 
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^—^—  ~      rviLiZ^j  j\.4.L2ZJ2  v  2     1  4     2  7     3/3  7~$ 

*    50j 


(K2Lj2  +   K4L22  +  K7L32)0!  +   K7L3L126; 


au 
*  ae2 


—  K.^L>j2^'3  KjLj2^4     "■      IVyljjL]^  j     "•      1V7LJ2   D 


2 


+   KALAL13cos(6o)02  ■    KALAL13cos(e0)e3 

au 
aea 


*  T7T=   "KALAL13cos(0o)e2  +   K^L^cos^o)©, 


'i 


4.        Generalized    Power    Lost    Due    to    Damping    Forces    (F) 

The    amount     of     generalized     power     that     is     dissipated     by 
each    of  the    dampers   in    the    system    given    by    the   equation: 

where: 

Fj  =    the    generalized    power    dissipated    in    each    damper 

B;  =    the    damping    constant    for    each    damper 

Ax;    =   the   difference   in   velocity   of  the  attaching  points   of  each 

damper. 

The  total  amount  of  generalized  power  dissipated  in  the  system  is 
found  by  summing  the  energy  dissipated  in  each  dampers.  Summing 
the    energy    dissipated    by    the    damping    forces    yields    the    following: 

F-^fc  -   Zj2  +  iB2(r2  -   Z2f  +  iB3(r4  -  Z3)2  +  \B<{^  -  02)2 

In     Lagrange's     equation     this     total     dissipated     energy,     F,     must      be 

differentiated   with   respect   to   the   velocities   of  each   of  the  generalized 

aF 
coordinates,     ■— -. 

aqi 
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Carrying    out    this    differentiation     and     again     linearizing     results     in     the 
following: 


*  aT  =  ^  "  B,Z3  +  BlLl01 

dF 

*  Jg~  "       B2Z2    "      B2Z3    "     B2L201 

*  dF=    _BlZl  "   B^2   +  <Bi+B2  +  ^)Z3  -   BjZ,  +  (-B.L,  +   B2L2)0, 

ap        •         • 

*  ttt  -   B3Z4  -   B3Z3  -  B3L30j  -  B3L1202 

dF 

— r-  -    Bjl^Zj  -  B2L2Z2  +  (B1L1  +   B2L2  +   B3L3)Z3  -  B3L3 

+  (B^2  +   B2L22  +   B^2)^  +   B.L.L^ 

dF  •  ■  2  ■  ■ 

*  T7T  "     B3L12Z3     ~     B3L12Z4    +     ^LsLl^l    +     (B»L12      +      B4)02     ~   B403 

•g-  -B4e2  +  B4e, 


5.       External    Forces    (Qi) 

The   external    forces    acting    on    the   system     are     the     roadway 
inputs      that     are     transmitted     through     the     tire     springs.  For     the 

computer  model  these  input  forces  are  treated  as  purely  vertical 
forces.  The  forces  result  from  the  road  inputs  causing  a  change  in 
length  of  the  tire  springs.  These  road  inputs  are  depicted  on  Figure 
2-10  and   defined   as   follows   : 


*  Ui  =   the   road   input   for   the   front   HUMMV   wheel 

*  U2  =   the  road   input   for   the   rear  HUMMV   wheel 

*  U3  =   the   road   input   for   the   howitzer   wheel. 
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6.         Assembling     Lagrange's     Equation 

Combining    the    parts    of    Lagrange's    equations    that    were 
developed    in    the    above    paragraphs    results    in    a    system    of   equations 
of  motion   in   the   form: 

[M]{x}  +  [C(x)]{x}  +  [K]{x}    =  [f] 

where: 

[M]   =  the  7x7  system   mass  matrix 

[C(x)]   =  the  7  x   7   system  damping   matrix   that  is  a  function  of 

whether   the   shock   absorbers   are   in   tension   or   compression 
[K]  =  the  7  x   7     system   stiffness  matrix 
{x},  {x},  &  {x}    =   7   x    1    column   vectors   that  contain   accelerations, 

velocities  and  positions  of  all   coordinates,   respectively, 
[f]    =   the   7x1    forcing   vector   that   results   from    the   road   inputs. 

The   [M],   [C(x)],    and    [K]    matrices    for     this     system     are     contained     in 
Appendix    B.         This   completes   the   development   of  the   heave   and 
pitch     EOMs.        The     resulting     EOMs     were     numerically    integrated      to 
determine    the    time    history    for    motions    of   the    component     masses    of 
the     system. 

C.       HUMMV/HOWITZER    MODEL    NUMERICAL    INTEGRATION 

The    following    section    describes    the    procedures    that    were    used 
to    numerically    integrate    the    previously    developed    equations    of 
motion. 

1.        Transformation     to     1st     Order     Equations 

In    order    to    numerically    solve    the    equations    of   motion,    the 
formulated    system    of    second    order    ordinary    differential    equations 
must    be    transformed    into    a    system    of    first    order    differential 
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equations.  This  equation  transformation  was  performed  using  matrix 
algebra  and  results  in  the  following  expression  for  the  system  of  first 
order     differential     equations: 

0 


,    ,   r  [o]       in    ifxi 

xdotj   =i  ,  I    .     + 

|_[-M~  K]    [-MMC(x)]  J  LxJ      [_M    f. 


where: 


[xdot]    =  a   14   x    1    column   vector  of  which   the   first   seven  elements 
are      the   velocities   of  all    the   generalized   coordinates   and    the    last 
seven    elements    are    the    accelerations    of    the    generalized    coordinates. 


=  a    14   xl    column   vector  of  which   the   first   seven   elements 


are    the   positions    in    terms    of   the    generalized    coordinates   and    the    last 
seven    elements    are    the    velocities    of    the    generalized    coordinates. 

[0]  =  a  7  x  7  zero  matrix 

[I]  =  a  7  x  7  identity  matrix 

[-M_1C(x)j  &[-M_1k]    =   these  are   7   x   7   matrices   that  result   from 

multiplying    the    negative    inverse    of   the    system    mass    matrix,    [M], 
times   the   damping   matrix,    [C],   and   the   stiffness   matrix,    [K], 
respectively 


=  a   14  x    1   vector  which  the    first    seven    elements    are    zero    and 


"0      " 

_M_1f_ 

the    last    seven    elements    are    the    matrix    product     of    multiplying     the     inverse 
of  the   mass   matrix,   [M],   times   the   forcing  vector,   [f]. 

For  the  integration  the  system  was  started  from  rest,  thus,  all 
seven  of  the  system's  displacement,  {x},  and  velocity,  {x},  state  variables 
had   a   value   of  zero   for  initial   conditions. 
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2.  Integration     Difficulties 

The  equations  of  motion  were  numerically  integrated  using  the 
Matlab  ODE  45  integration  routine.  This  is  a  variable  time  step  4th  order 
Runge-Kutta  integration  routine.  Integrating  these  equations  of  motion 
using  ODE  45  was  complicated  by  the  non-linear  damping  of  the  shock 
absorbers  and  the  required  shutoff  of  the  tire  springs  when  they  leave 
contact    with    the   road    surface. 

When  the  damping  values  of  the  shock  absorbers  would  abruptly 
change  from  compression  to  tension  and  vice  versa,  there  would  be  a 
corresponding  abrupt  change  in  the  shock  absorber  force.  Determining 
when  the  shock  absorber  damping  value  should  change  is  built  into  the 
integration  code  by  checking  for  each  time  step  if  each  of  the  shock 
absorber  forces  are  compressive  or  tensile.  An  "if"  statement  ensures  that 
the  correct  damping  value  is  used  for  the  condition  that  exists.  The  abrupt 
changes  in  force  in  the  shock  absorbers  slowed  down  the  tolerance  based 
ODE  45  routine.  A  similar  situation  occurred  when  the  tire  spring  forces 
would  abruptly  change  to  the  respective  static  weights  when  contact  with 
the   road   was   lost   for  a  particular   wheel. 

To  increase  the  speed  of  the  ODE  45  routine  the  integration  was 
divided  into  one  second  time  segments  so  the  solution  matrix  that  contains 
the     state     variables     remains     relatively     small.  The     displacements     and 

velocity  states  at  the  end  of  the  previous  time  segment  are  used  as  the 
initial    conditions    for    the    next   one    second    time    segment. 

The  resulting  time  histories  for  the  displacements  and  velocities 
allowed     post-processing     to     be     completed     determine     accelerations,     forces, 
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and    stress    in    key    system    components.        The    computer    code    developed    for 
this    numerical    integration    is    contained    in    Appendix    C. 
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III.       ROADWAY    INPUTS 

This  chapter  explains  the  road  inputs  that  were  used  to  drive 
the  computer  model.  The  roadway  inputs  are  defined  in  terms  of  the 
vertical  change  upward  or  downward  of  the  values  of  the  road  inputs 
Ui,  U2,  and  U3  that  are  shown  on  Figure  2-10.  These  are  independent 
vertical  displacement  inputs  to  tires  for  each  of  the  wheels  in  the 
computer  model.  As  defined  in  Chapter  II,  there  is  one  input  for  the 
front  HUMMV  wheels,  (Ui),  a  second  input  for  the  rear  HUMMV 
wheels,  (U2),  and  a  third  input  for  the  howitzer  wheels,  (U3).  The 
model  can  be  driven  over  any  terrain  that  can  be  defined  in  terms  of 
the  road  inputs  Ui,  U2,  and  U3  The  road  inputs  are  functions  of  time 
which  allows  the  user  to  control  the  speed  at  which  the  model  is 
driven  by  using  the  time,  velocity,  and  distance  relationships.  For 
this  design  study  the  road  inputs  used  were  a  washboard  course  that 
was  like  one  of  the  courses  used  in  the  actual  test  at  Aberdeen 
Proving  Grounds,  MD  (see  Chapter  VI),  and  a  predicted  "worst  case" 
pothole    road. 

A.      WASHBOARD  COURSE 

The  washboard  course  that  was  used  to  drive  the  model  was 
constructed  to  be  similar  to  the  Six  Inch  Washboard  Course  at  APG. 
During  the  instrumented  test  at  APG  there  were  several  data  runs 
taken  on  the  Six  Inch  Washboard  Course  (again,  see  Chapter  VI). 
Comparing    the    model   results    for   this    type   of   washboard    course   to    the 
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actual  instrumented  data  from  the  APG  test  for  a  similar  washboard 
provided  a  means  of  validating  the  output  results  of  the  computer 
model. 

To  model  the  APG  Six  Inch  Washboard  Course  a  sinusoidal  input 
was  used  with  an  amplitude  of  three  inches  and  a  distance  from  peak 
to  peak  of  72  inches.  However,  because  the  distance  between  the 
front  HUMMV  wheels  and  both  the  rear  HUMMV  wheels  and  also  the 
howitzer  wheels  does  not  correspond  to  even  multiples  of  the  72  inch 
peak  to  peak  distance,  a  lag  had  to  be  built  into  the  input  sinusoid  for 
both  the  rear  HUMMV  wheels  and  the  howitzer  wheels.  This  lag  was 
necessary,  to  get  the  correct  timing  of  the  bump  wheel  inputs  for  all 
the    wheels    in    the    system. 

Additionally,  in  order  to  prevent  difficulties  in  the  integration 
routine  that  would  occur  if  the  wheel  inputs  were  all  turned  on  at 
the  start  of  the  integration  run,  the  rear  HUMMV  and  howitzer  wheel 
inputs  were  not  turned  on  initially.  These  inputs  were  turned  on  at 
the  time  when  the  respective  wheels  would  start  hitting  the 
washboard  course  bumps  if  the  HUMMV/howitzer  system  entered 
the  course  from  level  terrain.  In  other  words,  the  front  HUMMV 
wheels  entered  the  washboard  and  started  taking  road  input,  Ui, 
while  the  rear  HUMMV  wheels  and  the  howitzer  wheels  remained  on 
level  terrain  which  provided  zero  input  to  them.  Then  the  rear 
HUMMV  wheels  entered  the  washboard  and  started  taking  road 
input,  U2,  while  the  front  HUMMV  wheels  continued  to  have  the 
input,     Ui,  and    the    howitzer    wheels  continued     to     have     zero     input. 
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Finally,    the    howitzer    wheels  entered    the    washboard    and    at    this    point 
all    three    wheels    had    their    separate    inputs    Ui,  U2,  and  U3. 

This  sinusoidal  input,  for  which  the  wavelength  was  defined  in 
terms  of  horizontal  distance  between  peaks,  was  converted  into  a 
function  of  time  by  using  time,  velocity,  and  distance  relationships. 
To  calculate  the  time  between  successive  peaks  for  a  specified 
velocity  the  inputs  Ui,  U2,  and  U3  were  then  defined  in  terms  of  time 
which  allowed  the  computer  model  to  be  "driven"  at  any  desired 
speed.  Making  the  input  a  function  of  time  was  required  for  using 
the  time  based  ODE  45  integration  routine  of  Matlab  to  solve  the 
equations  of  motion.  Figure  3-1  depicts  the  wheel  inputs  Ui(  U2,  and 
U3  for   the   washboard    course   being    driven    at    5.0   MPH. 
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Solid    Howitzer  Wheel  Input,  U3 


Figure    3-1       Six    Inch    Washboard   Course    Wheel    Inputs 
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B.     "WORST  CASE"  POTHOLE  COURSE 

The  model  was  also  "driven"  over  a  trench  type  pothole  to 
simulate  some  of  the  most  severe  road  conditions  that  the  howitzer 
would  be  subjected  to  during  tactical  operations.  The  model  results 
for  A-Frame  force  from  this  terrain  were  used  to  predict  if  the  A- 
Frame  and  T-Bar  were  strong  enough  to  handle  this  "worst  case" 
terrain.  The     pothole     course     road     inputs,     Ui,    U2,     and    U3,    are 

determined  by  a  separate  computer  program  (see  Appendix  D)  that 
computes  the  changes  in  axle  height  as  each  of  the  HUMMV  wheels 
and  the  howitzer  wheels  roll  through  the  pothole.  The  road  inputs 
are  based  on  the  path  that  would  be  traced  out  by  the  center  of  a  17 
inch     radius     wheel      rolling     through     the      specified      pothole.  The 

computer  program  that  was  developed  to  determine  the  road  inputs 
allows  the  user  to  specify  the  desired  dimension  of  the  pothole  for 
driving  the  computer  model.  For  this  design  study  the  "worst  case" 
pothole  was  determined  to  be  20  inches  deep  and  30  inches  across 
with  85  degree  slopes  on  each  edge.  This  size  was  selected  because  it 
was  judged  to  be  the  largest  pothole  that  the  HUMMV  driver  would 
traverse  without  stopping  the  HUMMV.  Thus,  the  "worst  case" 
pothole  dimensions  were  selected  based  on  what  was  judged  to  b  e 
the  most  severe  terrain  the  HUMMV/howitzer  would  encounter 
while  traveling  during  tactical  operations.  The  "worst  case"  pothole 
selected  was  slightly  larger  than  the  dimensions  of  the  pothole  that 
was  used  in  an  informal  test  conducted  at  RIA  in  June  1992.  Figure 
3-2  shows  the  pothole  profile,  and  the  axle  path  traced  for  the 
howitzer   or    HUMMV    wheel    rolling    through    this    "worst   case"    situation. 
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Figure  3-2      "Worst  Case"   Pothole  Profile  and  Axle  Trace 
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IV.      MODEL   RESULTS 

This  chapter  contains  the  results  of  driving  the  computer  model 
over  the  terrain  that  was  explained  in  Chapter  III.  The  results  from 
driving  the  model  over  the  washboard  course  show  close  agreement 
actual  results  obtained  from  the  APG  test.  A  detailed  comparison  of 
the  model  results  to  the  actual  results  is  contained  in  Chapter  VI. 
The  model's  results  from  the  "worst  case"  pothole  were  used  to 
predict  the  maximum  value  of  A-Frame  force  for  the  howitzer.  From 
this  predicted  maximum  force  the  stress  in  the  A-Frame  and  T-Bar 
was  calculated  and  compared  to  the  yield  stress  for  the  material  (see 
Chapter    V). 

The  predicted  A-Frame  force  is  calculated  by  multiplying  the 
combined  spring  stiffness  of  the  A-Frame/T-bar  times  the  change  in 
length  of  A-Frame.  The  change  in  length  of  the  A-Frame  is  a 
function  of  the  difference  in  the  angular  displacements  of  the  gun 
tube  mass  and  the  howitzer  trail/saddle  mass  at  a  given  time.  The 
angular  displacements  of  the  gun  tube  and  the  trail  were  determined 
from    numerically    integrating    the    equations    of    motion. 

In  addition  to  determining  results  in  terms  of  A-Frame  force, 
the  system's  mass  and  stiffness  matrices  (see  Appendix  B)  were  used 
to  determine  natural  frequencies  and  mode  shapes  of  the  system. 
These  natural  frequencies  and  mode  shapes  were  used  to  gain 
physical      insight      into     the     system's      performance.  The      natural 

frequencies    and    mode    shapes    for    the    system    are    shown    in    Appendix 
E 
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The  Matlab  "FFT"  function  was  used  to  perform  the  frequency 
analysis  on  the  A-Frame  force  results  from  the  washboard  course  to 
determine  the  frequencies  at  which  the  A-Frame  force  occurred  in 
the  model  results.  This  frequency  analysis  a  used  1250  Hz  data 
sampling  rate.  This  sampling  rate  yielded  a  time  increment,  dt,  that 
equals  1/1250  or  the  inverse  of  the  sampling  rate.  In  order  to  have 
a  frequency  increment  that  was  sufficiently  small,  8192  data  points 
were  used  for  the  "FFT".  The  record  length,  T,  equals  dt  times  the 
number  of  points  or  8192/1250.  The  frequency  increment,  df,  was 
the  inverse  of  the  record  length  or  1250/8192. 
A.      WASHBOARD  COURSE  RESULTS 

The  time  histories  of  the  computer  model's  predictions  for  A- 
Frame  force  from  the  Six  Inch  Washboard  Course  are  periodic  with  a 
period  that  corresponds  to  the  period  of  the  washboard  input.  But 
unlike  the  sinusoidal  washboard  input  the  A-Frame  force  output 
predictions  are  not  sinusoidal.  These  non-sinusoidal  A-Frame  force 
outputs  most  likely  results  from  the  effects  of  the  non-linear  shock 
absorber    damping    that    was    incorporated    into    the    model    (see    Chapter 

II). 

The  frequency  analysis  of  the  computer  model's  predictions  for 
A-Frame  force  shows  that  as  expected  the  largest  component  of 
output  force  occurs  at  the  same  frequency  as  the  washboard  input 
frequency.  In  all  cases,  there  is  also  a  significant  but  lesser  force 
component  that  occurs  at  a  frequency  that  is  twice  the  input 
frequency. 

The    computer    model    was    driven    over    the     six     inch     washboard 

course    at    speeds    of    3.75,    5.0    MPH,    and    5.25    MPH.       These      speeds 
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correspond  to  the  speeds  for  which  actual  data  was  taken  during  the 
APG  test.  The  A-Frame  force  levels  were  relatively  low  at  3.75  MPH. 
The  force  levels  in  the  A-Frame  reached  a  maximum  at  5.25  MPH 
because  this  speed  for  the  washboard  course  corresponds  to  a 
natural  frequency  of  the  system.  Driving  the  model  at  speeds  higher 
than  5.25  MPH  resulted  in  decreasing  A-Frame  force  levels.  The 
time  histories  and  frequency  analysis  of  the  model's  predicted 
results    for    A-Frame    force   are    shown    in    Figures    4-1    through    4-6. 

1.  Results    at    3.75    MPH 

The  time  domain  model  results  at  3.75  MPH  are  shown  in 
Figure  4-1.  The  predicted  A-Frame  forces  are  approximately  1500 
lbs  in  tension  and  1200  lbs  in  compression.  Figure  4-2  shows  the 
frequency  analysis  at  this  speed.  As  expected,  there  is  a  large  force 
component  at  the  input  of  3.75  MPH  or  a  frequency  of  0.92  Hz  and 
there   is    smaller   force   spike   at   twice   this   frequency,    1.84   Hz. 

2.  Results    at    5.0    MPH 

The  time  domain  results  are  shown  in  Figure  4-3.  The  A- 
Frame  force  levels  reach  2800  lbs  in  tension  and  4800  lbs  in 
compression.  Figure  4-4  shows  the  frequency  analysis  results.  Again, 
a  large  force  component  at  the  input,  5.0  Hz,  and  another  force  spike 
at   twice   this    frequency,    2.42   Hz. 

3.  Results    at    5.25    MPH 

The    time    domain    results    are    shown    in    Figure   4-5.      The   A- 

Frame     force     levels     reach     2900     lbs     in     tension     and     4700     lbs     in 

compression.        Figure    4-6    shows    the     frequency     analysis    results.        A 

large    force    component    exists    at    the    input    frequency,     1.28    Hz,     and 

smaller    force    component   exists    at    twice    this    frequency,    2.56    Hz. 
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Figure    4-1        Model    Results    for    A-Frame    Force    at    3.75     MPH     (Time 
History) 
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Figure  4-2      Model   Results   for  A-Frame  Force  at  3.75   MPH   (FFT) 
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A -Frame  Force  5.0  MPH 
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Figure    4-3        Model     Results     for    A-Frame    Force    at    5.0    MPH     (Time 
History) 
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Figure  4-4      Model   Results   for  A-Frame  Force  at  5.0  MPH   (FFT) 
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Figure    4-5        Model    Results    for    A-Frame    Force    at    5.25     MPH     (Time 
History) 
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Figure  4-6     Model   Results   for  A-Frame  Force  at  5.25   MPH   (FFT) 


49 


B.  NATURAL   FREQUENCIES   AND   MODE   SHAPES 

The  natural,  frequencies  and  mode  shapes  of  the  system  were 
determined  by  using  the  mass  and  stiffness  matrices  that  were 
developed  in  Chapter  II.  The  eigenvalues  and  eigenvectors  of  the 
system  are  determined  from  the  matrix  that  results  from  multiplying 
the  inverse  of  the  mass  matrix,  [M-  *],  by  the  stiffness  matrix,  [K]. 
The  system's  natural  frequencies  are  the  square  root  of  the 
eigenvalues.  The  system's  mode  shapes  are  the  eigenvectors  that 
correspond  to  these  natural  frequencies.  Analysis  of  these  natural 
frequencies  and  mode  shapes  shows  that  the  howitzer's  bounce  mode 
occurs  at  a  natural  frequency  of  1.28  Hertz  which  corresponds  to  a 
speed  of  5.25  MPH  for  the  howitzer  traveling  over  the  washboard. 
By  providing  an  excitation  at  a  natural  frequency  of  the  system  (1.28 
Hz)  the  maximum  A-Frame  force  that  could  possibly  be  developed 
for    this    particular    washboard    input    resulted. 

C.  "WORST  CASE"  POTHOLE  COURSE  RESULTS 

The     model     results     from     the     "worst     case"     pothole     that      was 

identified    in    Chapter    III    show     that     maximum     A-Frame     force    occurs 

lbs  at   15  MPH.     At  speeds  slower  or    faster    than    this    speed    the    model 

results    showed   lower   levels   of  A-Frame   force   for   this   pothole.      Figure 

4-7    shows   the   A-Frame   force   for   the   "worst   case"   pothole   at    15    MPH. 

There    are     spikes     in     the     A-Frame     force     when     the    wheels     of     the 

system    hit    the    pothole.        On    Figure    4-7,    at    0.1     seconds    the    front 

HUMMV  wheels    hit,    at    0.6    seconds    the    rear   HUMMV    wheels    hit,    and 

at     1.2     seconds     the     howitzer    wheels     hit     the     pothole.         When     the 
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howitzer   wheels   hit   at    15    MPH    the   maximum     A-Frame     force     is     8,500 
lbs. 
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Figure  4-7      Model   Results    of   A-Frame    Force    for    "Worst    Case"    Pothole 
at    15   MPH 
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V.       STRUCTURAL    ANALYSIS 

This  chapter  explains  the  structural  analysis  that  was 
performed  on  the  howitzer  in  order  to  determine  both  the  force 
levels  in  the  A-Frame  and  the  stress  on  the  howitzer's  components. 
To  determine  the  actual  force  levels  in  the  A-Frame,  the  A-Frame 
was  instrumented  with  strain  gages  as  shown  in  Figure  5-1.  This 
instrumented  A-Frame  was  then  calibrated  to  determine  a 
calibration  factor  which  would  correlate  strain  gage  readings  to  A- 
Frame  force  levels.  This  calibrated  A-Frame  was  used  during  the 
dynamic  test  (see  Chapter  VI)  at  Aberdeen  Proving  Grounds  (APG)  to 
measure    the    actual    force    levels    in    the    A-Frame. 

A.       A-FRAME    CALIBRATION 

The  primary  purpose  of  calibrating  the  A-Frame  was  to  develop 
a  reliable  means  to  correlate  dynamic  strain  gage  measurements  on 
the  legs  to  the  amount  of  force  that  causes  the  strain.  Additionally, 
the  strain  measurements  obtained  during  calibration  were  used  to 
determine  the  value  of  the  equivalent  linear  spring  stiffness  of  the 
A-Frame  for  the  computer  model.  During  the  calibration  extra  strain 
gages  were  attached  to  help  identify  levels  of  stress  at  key  locations 
other  than  the  locations  needed  for  the  force  versus  strain 
calibration.  These  extra  gages  provided  information  that  was  used  to 
form    conclusions    about    the    stress    in    the    A-Frame,    discussed    later. 
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ROSETTE 
G11.G12 
&G13 


NOTE:  There  were  ten  single  strain  gages,  Gl  through  G10, 
located  as  indicated  on  the  diagram.  For  the  rosette,  Gl  1 
is  the  transverse  leg,  Gl  2  is  the  45  degree  leg  and  Gl  3  is 
the  axial  leg. 


Figure   5-1       Location   of  A-Frame   Strain    Gages 
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The  A-Frame  was  loaded  in  both  a  vertical  and  lateral  load 
configuration  to  determine  if  there  were  any  significantly  different 
effects  between  the  two  loading  configurations.  Figures  5-2  &  5-3 
show    the    test    setup    for    both    these    configurations. 

1.  A-Frame     Instrumentation 

To  conduct  the  calibration,  the  A-Frame  was  instrumented 
with  13  strain  gages  (Figure  5-1).  The  data  taken  from  the  pair  of 
single  strain  gages  located  on  the  bottom  of  each  leg  was  used  for  the 
development  of  the  calibration  curve.  The  pair  of  single  strain  gages 
on  the  top  of  each  leg  were  used  to  verify  that  loads  in  the  legs  were 
mainly  axial  loads  as  opposed  to  bending  loads.  The  single  gage  on 
the  A-frame  cross  support  brace  and  the  single  gage  on  the  attaching 
clamp  for  this  brace  were  used  to  identify  stress  levels  in  these 
members.  The  latter  two  strain  gages  were  used  because  there  were 
some  known  cases  where  the  cross  support  brace  clamps  had 
cracked,  and  it  was  unknown  whether  these  cracks  were  the  result  of 
high    load    levels    or    caused    by    manufacturing    defects. 

2.  Support     Equipment 

The  calibration  was  conducted  using  the  lab  facilities  of  the 
Department  of  Mechanical  Engineering  at  the  Naval  Postgraduate 
School.  To  hold  the  A-Frame  during  the  calibration  a  test  support 
structure  was  built.  The  main  objective  of  this  support  stand  was  to 
closely  match  the  mechanical  effects  that  occur  in  the  A-Frame  when 
it  is  loaded  while  connected  to  the  T-Bar.  The  test  support  used  was 
constructed  out  of  heavy  steel  and  used  actual  eyebolts  from  a  Ml  19 
Howitzer   to   connect   the   bottom   of  the   A-Frame   legs   to   the   test 
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Figure    5-2       Vertical    Configuration    for    A-Frame    Test 


Figure    5-3       Lateral    Configuration    for    A-Frame    Test 


5  5 


support  (Figure  5-4).  A  clamping  device  to  attach  the  apex  of  the  A- 
Frame  to  the  load  cell  was  also  designed  which  simulated  the  loading 
effects  of  the  actual  clamping  device  on  the  howitzer  (Figure  5-5). 
The  test  support  structure  achieved  the  goal  of  simulating  the  A- 
Frame   being    attached   to    the   T-Bar   of   the    howitzer. 

3.  Load    Cells 

The  vertical  loads  were  applied  using  a  Riehe  Testing 
machine  as  shown  in  Figure  5-2.  This  loading  device  was  used  to 
apply  vertical  loads  up  to  12,000  lbs.  The  lateral  loads  were  applied 
using  a  turn-buckle  that  was  connected  to  a  5,000  pound  rated 
capacity  dynamometer.  This  dynamometer  and  turn-buckle  system 
is   shown   in   Figure   5-3. 

4.  Calibration     Procedures 

Two  data  runs  were  recorded  for  both  the  vertical  and 
lateral  configurations.  A  data  run  consisted  of  loading  the  A-Frame 
in  increments  to  the  maximum  load  for  that  particular  data  run  and 
then  unloading  it  in  the  same  manner.  All  strain  gages  were  zeroed 
before  the  start  of  each  run.  At  each  load  increment,  strain  gage 
readings  were  recorded  for  all  of  the  strain  gages.  Tables  F-l 
through  F-4  contain  the  strain  gage  readings  and  load  data  obtained 
from    these    calibration    procedures. 
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Figure    5-4       A-Frame    Connection    to    Test    Mount 


Figure    5-5       A-Frame    Apex    Clamping    Device 
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5.        Calibration     Factor 

From    the    vertical    test    strain    gage    data    in    Appendix    F  the 

A-Frame       vertical       force       versus       strain       calibration       factor  was 

determined.        Since    the    A-Frame    vertical    load    test    was    a    static  test 

the  vertical  force  loading  at  the  apex  of  the  A-Frame  must  be 
balanced  by  the  sum  of  the  vertical  reaction  forces  in  the  legs  of  the 
A-Frame       (Figure    5-6). 


\ 


Strain 
Gage  Pair    #1 


in 

Pair  #  4 


where: 

F  =  the  applied    vertical  load  at  the  apex  of  the  A-frame 

Rhi  &  ^2   =  tne  horizontal  reaction  forces  where  the  A-frame 
connects  to  the  T-BAR 

Rvl  &  1^2    =  tne  vertical  reaction  forces  where  the  A-frame 
connects  to  the  T-BAR 


Figure   5-6      Balance  of  Forces  on   the   A-Frame 
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The  reaction  forces  (Rv  1  &RV  2)  at  the  bottom  of  each  leg  of 
the  A-Frame  were  determined  by  using  the  strain  gage  readings 
from  Gage  Pair  1  (gages  1  &  2)  and  Gage  Pair  4  (gages  9  &  10)  that 
are  identified  in  Figure  5-6.  The  average  value  of  the  strain  for 
either  of  the  strain  gage  pairs  was  multiplied  by  the  Modulus  of 
Elasticity  of  the  A-Frame  (10  x  106  psi)  which  yields  the  axial  stress 
in  the  A-Frame  leg.  This  axial  stress  value  was  then  multiplied  by 
the  cross  sectional  area  of  the  A-Frame  legs  (0.6107  in2)  to 
determine  the  axial  force.  To  calculate  the  vertical  component  of  this 
axial  force  the  above  value  for  axial  force  was  multiplied  times  the 
cosine  of  28°  A  sample  calculation  is  shown  below  using  the  10,000 
lb  load  increment  of  Vertical  Test  #1  (see  Table  F-l  for  strain  gage 
data): 


Strain  Gage  Pair  1    =>  Avg   Strain  -  988  *  10"6 


.-.  Stress  =  Strain  x  Modulus  of  Elasticity 
=  (988  x  10"6)(l0  x  106psi)  =  9,880  psi 
=»  Force  =  Stress  x  Area  =  (9,880  psi)(0.6107  in2)  =  6,033  lbs 
.*.  The  vertical    reaction  force:   Rvl  =  (6,033  lbs)(cos(28))  =  5,327  lbs 

Strain  Gage  Pair  4    =>  Avg  Strain  =  1,095  *  W6 
.-.  Stress  -  Strain  x  Modulus  of  Elasticity 

=  (1,095  *  10"6)(l0  *  106 Psi)  =   1,095  psi 

=>  Force  =  Stress  *  Area  =  (10,950  psi)(0.6107  in2)  =  6,687  lbs 

.-.The  vertical    reaction  force:   Rv2  =  (6,687  lbs Xcos( 28))  =  5,904  lbs 

.£  Vertical   Forces   =►  Rvl  +  Rv2  =  5,327  +  5,904  =  11,231  lbs 
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To  determine  the  calibration  factor  for  vertical  loading, 
the  percentage  difference  between  the  known  apex  load  and  the  sum 
of  the  calculated  vertical  force  components  at  the  bottom  of  the  A- 
Frame  was  determined.  The  average  percent  difference  for  all  the 
load  increments  is  equal  to  the  calibration  factor.  This  calibration 
factor  can  then  be  used  to  compare  measured  loads  in  the  A-Frame 
legs  to  the  applied  load  at  the  apex  of  the  A-Frame.  An  average 
calibration  factor  of  0.88  was  determined  over  all  of  the  load 
increments    from     both    vertical    tests.  A    sample    calculation    of    the 

calibration  factor  for  the  10,000  lb  force  increment  of  Vertical  Test 
#1    (see  Table  F-l)   is   shown   below: 

Actual  Vertical   Force   =>  10,000  lbs 
Measured  Vertical   Force    =>  11,231  lbs 

Calibration   Factor   =>      ,K^  =  0.89 

11,231 

6.         A-Frame     Stiffness 

The  measured  strain  data  was  also  used  to  determine  a  n 
effective    linear    spring    stiffness    for    the    A-Frame.  To    determine    this 

stiffness  the  applied  vertical  load  for  a  specified  load  increment  was 
divided  by  the  corresponding  amount  of  vertical  deflection  in  the  A- 
Frame.  The  average  axial  deflection  in  the  A-Frame  legs  was 
calculated  from  the  average  strain  at  the  bottom  of  each  leg.  The 
average  axial  deflection  was  then  multiplied  by  the  cosine  of  28°  to 
calculate  the  amount  of  vertical  deflection.  A  sample  calculation  for 
the  10,000  lb  load  increment  of  Vertical  Test  #1  (see  Table  F-l)  is 
shown    below: 
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Strain  Gage  Pair  1  =  998;       Strain  Gage  Pair  4  =  1095 

Avg   axial   strain  = =  1042  *  10"6;    Multiply  by  leg  length 

=>  Avg  axial   deflection  =  (28  in)(l042  *  10"*)  =  0.029  in 
Convert   to  vertical    deflection    =►  (0.029  in  Xcos(28))  =  0.0257  in 

.-.  A  -  Frame  Stiffness  =  10'000  lbs  =  388,400  lbs/  in 

a  0257  in 


7.         Observations 

From    analyzing    the    results    of    both    the    vertical    and     lateral 
A-Frame    loadings    the    following    conclusions    were    reached    regarding 
A-Frame     loading: 

*  The  strain  measurements  in  both  the  cross  member  brace  and 
the  attaching  clamp  for  this  brace  on  the  A-Frame  (strain  gages  3 
&  4,  in  Figure  5-1)  were  very  low  regardless  of  the  force  at  which 
the  A-Frame  was  loaded.  This  indicates  that  these  members  are 
not    critical    components   of  the   A-Frame    in    regards    to    failure. 

*  The  A-Frame  loads  are  carried  mostly  through  axial  loads  in  the 
A-Frame  legs  as  opposed  to  bending  loads.  This  means  that  the 
bottom  connection  joints  where  the  A-Frame  attaches  to  the  T-Bar 
act  as  pinned  joints  However,  because  there  was  some  bending  in 
the  A-frame  legs  at  lower  force  levels  the  bottom  joints  do  not  act 
as  perfect  pin  joints.  The  higher  force  levels  break  the  static 
friction  effects  in  the  joints  causing  the  joints  to  act  more  as 
perfect    pin    joints. 

B.       STRESS    ANALYSIS 

The    stress    analysis     focused     on     determining     the     highest     stress 
levels    in    the    T-Bar    and    the    A-Frame.  To    determine    if    mechanical 

failure    will    occur    due    to    exceeding    the    yield    stress,    the    stress   at     the 
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probable     points     of    highest     stress     in     the     A-Frame     and     T-Bar    were 
compared    to    the    yield    stress    values    for    the    respective    materials. 
1.        T-Bar    Stress 

The  T-Bar  and  A-Frame  assembly  are  shown  below  in 
Figure  5-7.  Point  A  on  this  figure  is  identified  as  the  point  of  highest 
stress.  This  is  where  the  geometry  of  the  T-Bar  assembly  changes 
and  the  stronger  built  up  portion  of  the  T-Bar  ends.  To  simplify  the 
stress  analysis  the  T-Bar  is  modeled  as  a  cantilevered  beam.  Figure 
5-8  below  depicts  the  cantilevered  beam  modeling  of  the  T-Bar  as 
viewed  from  the  end  of  the  T-Bar  and  the  rear  of  the  howitzer.  The 
point  of  highest  stress  was  determined  to  be  at  the  corner  of  the  T- 
Bar  at  the  base  of  the  modeled  cantilevered  beam.  One  of  these 
corners   is   identified   as   Point   A   in   Figure   5-8. 


A-Frame 


Transverse 
Hole 


High  Stress 
Location 
(Point  A) 


High  Stress 
Location 
(Point  A) 


8  inches 


8  inches 


^T-B 


ar 


Figure    5-7      T-Bar    and    A-Frame    Assembly 
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View  from  Rear  of  Howitzer 


View  from  end  of  T-Bar 

z 

p        i 


Point  A 


Figure   5-8         T-Bar   Cantilevered   Beam   Model 

To  calculate  the  stress  at  Point  A,  the  normal  stress 
components  that  result  when  the  axial  force  in  the  A-Frame  leg,  P, 
was  resolved  into  components  in  the  x,  y,  and  z  directions  as  shown 
on  Figure  5-8.  The  stress  at  Point  A  is  the  sum  of  these  .normal  stress 
components  using  the  principle  of  superposition  of  stress.  To  find 
the  normal  stress  components  first  the  A-Frame  leg  force,  P,  was 
projected  into  a  force  component  along  the  axis  of  the  T-Bar  in  the  y 
direction,  Py  =  sin(28)P,  and  a  second  force  component  in  the  xz 
plane,  Pxz=  cos(28)P.  The  Pxz  component  was  further  broken  up  into 
components  in  the  x  and  z  directions:  Px  =  cos(27)Pxz  and  Pz  = 
sin(27)Pxz.  To  calculate  the  normal  stress  that  was  caused  by  these 
force  components  Px,  Py,  and  Pz  well  known  strength  of  materials 
formulas  for  stress  in  a  beam  are  applied.  Therefore,  applying 
these  formulas  the  expression  for  the  maximum  stress  in  the  T-Bar 
is: 
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G  = 

where: 


'PXxN 


I 


(Eqn  5-1) 


Px,  Py,  and  Pz:      are   defined   in    the   above   paragraph. 

L=8   inches:      The    distance    from    the    base    of   the    cantilever    to    the 
point    of    application. 

A=.9240    in2:      The  cross  sectional   area  of  the  T-Bar. 

Ix=.8962    in4:      The  bending   moment   of  inertia  about   the   x   axis. 

Iz=.5979    in4:      The   bending   moment   of  inertia  about   the   z   axis. 

x   =1.300  in:      The  distance  in   the  x   direction   from   the  geometric 
center  of  the  T-Bar  to   the  outer  edge  of  the  T-Bar. 

z   =1.005   in:      The  distance  in   the  z  direction   from   the  geometric 
center  of  the  T-Bar   to   the  outer  edge   of  the  T-Bar. 


2.        A-Frame    Stress 

The  probable  point  of  highest  stress  in  the  A-Frame  will 
occur  around  the  0.375  inch  diameter  transverse  circular  hole  that 
positions  the  cross  member  attaching  clamps.  This  hole  is  identified 
in     Figure     5-7     above.  At     the     location     of     this     hole     a     stress 

concentration  will  exist.  Figure  5-9  below  contains  the  empirical 
formula  from  page  601  of  [Ref.  4]  that  was  used  to  determine  the 
stress    concentration    in    a    pipe    section    with    a    transverse    hole. 

From  the  above  formula  the  stress  concentration  factor  for 
axial      loads      was      determined      to      be      3.73.  Additional      stress 

concentrations  due  to  the  effects  of  bending  and  torsional  loads  were 
neglected.      The   justification    for   neglecting    the   effects   of  bending 
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fladiaJ  hole  in  a  hollow  or  solid      i  18a     FJasiK   iu»i,  axiaJ  icmion 
circular  snail 


- = >ij^-)  -he" ' = *' + *(S) + M£)s + *(£)' 


HD2 
j  and  where  for  HD  <  0.9  and  ir/D  <  0.45 

i  A,  =       3.000 

I  A.  =       2.773  +  1.529rf/D  -     4.379(d/£>)2 

A',  =   -0.421  -  \11SU/D  +  22.781(rf/D)J 

A,  =     16.841   +  16.678<//£>  -  4OOO7(rf/0)2 


6    D 

L 


For  a  solid  shafi  </  =  0 
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18b.  FJasuc  stress    bendinc  when  hole  o_      =  i : where  i  =  A",  +  A", I  —  I  +  A\( —  )    +  A. {  —  ) 

,     ■        ,  /  ^(O4  -  </4)  \x7/  VD/  '\D/ 

is  farthest  bom  bending  axis  '  ' 

and  where  for  HD  <  0.9  and  ir/D  <  0.3 

A,  =         3.000 

A,  =  -6.690-  1.62O<//0+  4.432(y/£>)2 
A",  =  44.739  +  I0.724<//D  -  19.927(</7£>)2 
A4  =  -53.307  -  2S.998J/D  +  43.258(<//D)2 


18c     Elastic  stress,  torsion 
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and  where  for  d/D  <  0.9  and  2r/£>  <  0.4 

A',  =  4.000 

A,  =  -6.793  +  1 A33J/D  -  0.126(rf/D)2 
A,  =  38.382  -  l.HUID  +  6.495<<//£>)2 
A4  =   -44.576  -  7  428*70  +  58.656<<i7D)2 


Figure    5-9       Stress    Concentration    Due    to    a    Circular    Hole    in    a    Pipe 
Section 


loads  was  that  both  the  A-Frame  calibration  results  (see  Appendix  F) 
and  the  results  of  the  instrumented  test  (see  Chapter  VI)  show  that 
the  joint  where  the  A-Frame  connects  to  the  T-Bar  acts  as  a  nearly 
perfect  pin  joint  at  higher  load  levels.  Any  stress  concentration 
reduction  due  to  the  spring  pin  that  is  driven  into  the  hole  were 
neglected  because  this  stress  reduction  was  assumed  to  be  small. 
Therefore,  maximum  total  stress  in  the  A-Frame  was  based  on  only 
the  axial  stress  component  in  the  legs  and  was  determined  using  the 
following    formula: 


a  = 


PK 


(Eqn  5-2) 


where: 

P  =     the  axial  load  in  the  A-Frame  leg. 

K  =  3.73:      the  stress  concentration  factor    for   axial    loading    that   is 
caused   by   the   hole   in   the   A-Frame   leg. 

A  =  0.6107  in2:      the  cross  sectional  area  of  the  A-Frame  leg. 
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Without  this  hole  in  the  A-Frame  the  maximum  allowable 
stress  would  be  increased  by  the  amount  of  the  stress  concentration 
factor,  3.73.  This  hole  appears  to  only  aid  in  assembling  the  A-Frame 
and  serves  no  other  function.  Elimination  of  this  hole  would  make 
the    factor    of   safety    more    than    three    times    greater. 

3.        Worst    Case    Stress    and    Factors    of    Safety 

Using  Eqn  5-1  and  Eqn  5-2  above,  the  maximum  stress  was 
calculated  for  the  T-Bar  and  the  A-Frame.  The  maximum  predicted 
A-Frame  force  of  8,500  lbs  from  driving  the  computer  over  the 
"worst  case"  pothole  described  in  Chapter  II  was  divided  by  two  to 
calculate  the  force  in  each  leg  of  the  A-Frame.  The  resulting  axial 
force,  P,  in  each  leg  of  the  A-Frame  equals  4,250  lbs.  From  Eqn  5-1 
the  resulting  stress  in  the  T-Bar  is  63,900  psi.  From  Eqn  5-2  the  A- 
Frame    stress    equals    26,000    psi. 

The  above  stress  levels  are  used  to  calculate  the  predicted 
factors  of  safety  for  the  T-Bar  and  the  A-Frame.  The  factor  of  safety 
is  calculated  by  dividing  the  yield  stress  of  the  component  by  the 
predicted  maximum  load  stress.  The  yield  stress  for  the  T-Bar  and 
A-Frame  are  116,000  psi  and  50,000  psi,  respectively.  Based  on  the 
maximum  predicted  stress  levels  given  above,  the  factor  of  safety  is 
1.81  for  the  T-Bar  and  1.92  for  the  A-Frame.  These  factors  of  safety 
show  that  both  the  T-Bar  and  the  A-Frame  are  strong  enough  to 
withstand    the    stress    from    the    worst    case    predicted    loads. 
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VI.       INSTRUMENTED    TEST 

An  instrumented  dynamic  test  was  conducted  on  the  Ml  19 
howitzer  at  Aberdeen  Proving  Grounds  (APG)  in  April  of  1993.  The 
howitzer  was  instrumented  as  described  below  in  paragraph  A.  The 
dynamic  testing  was  conducted  by  driving  the  howitzer  over  several 
of    the    APG    courses    at    various     speeds.  The     calibrated     A-Frame 

described  in  Chapter  V  was  used  to  measure  the  force  levels  in  the 
A-Frame  during  the  test.  In  most  respects,  the  test  results  for  A- 
Frame  force  are  similar  to  the  results  obtained  from  the  computer 
model.  However,  because  certain  physical  effects  were  not  included 
in  the  computer  model,  there  are  some  differences  between  the 
results  of  the  actual  test  and  the  results  from  the  computer  model. 
The  results  of  this  instrumented  test  confirmed  that  the  M 1  1 9 
Howitzer  could  be  safely  transported  in  the  firing  position  with  the 
A-Frame     attached. 

A.        INSTRUMENTATION 

For  this  test  the  howitzer  was  instrumented  with  20  strain 
gages,  two  tri-axial  accelerometers,  and  two  extensiometers.  This 
instrumentation    is    explained    in    the    following    paragraphs. 

1.  Accelerometers 

The  howitzer  was  instrumented  as  shown  in  Figure  6-1.  A 
tri-axial  accelerometer  was  placed  on  the  T-Bar  to  measure  the 
vertical,    longitudinal,    and    transverse    accelerations    of    the    howitzer 
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Figure    6-1       Howitzer    Instrumentation    for    APG    Test 
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trail  mass.  This  accelerometer  was  not  placed  at  the  center  of 
gravity  of  the  trail  mass  because  of  difficulties  in  attaching  an 
accelerometer  at  this  point.  The  location  chosen  on  the  T-Bar 
provided  accelerations  on  the  same  order  of  magnitude  and 
representative  of  the  accelerations  at  the  center  of  gravity  of  the 
trail.  The  other  tri-axial  accelerometer  was  attached  on  the  elevating 
mass   near   its  center  of  gravity.      This   location   is   shown   in   Figure   6-2. 

2.  Single    Strain    Gages 

The  two  pair  of  single  strain  gages  on  the  bottom  of  each  leg 
(gage  pair  1  &  4  in  Figure  5-6)  of  the  previously  calibrated  A-Frame 
were  used  to  measure  the  strain  in  the  A-Frame  legs.  The  calibration 
factor  determined  in  Chapter  V  was  used  to  correlate  measured 
strains  in  the  A-Frame  legs  to  force  levels  in  the  A-Frame.  Single 
strain  gages  were  also  attached  on  the  top  and  the  bottom  of  the  T- 
Bar  on  both  sides  of  the  howitzer.  The  location  of  these  gages  is 
shown  in  Figure  6-3.  Because  the  highest  levels  of  stress  on  the  T- 
Bar  are  at  the  corners,  it  was  not  possible  to  place  the  gages  at  these 
points.  The  gages  were  placed  as  close  to  the  high  stress  points  as 
possible.  Even  though  these  gages  could  not  be  used  to  determine 
the  maximum  stress  on  the  T-Bar,  they  could  be  used  to  roughly 
judge   the   level   of  stress. 

3.  Other    Gages 

Strain  rosettes  were  placed  on  the  bottom  of  the  A-Frame 
legs   and   on   the   bottom   howitzer   trail   on   each   side   of  the   howitzer. 
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Figure    6-2       Location    of   Elevating    Mass    Accelerometer 
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Figure   6-3      Location   of  T-Bar   Strain   Gages 
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An  extensiometer  was  used  on  each  shock  absorber  to  determine  if 
the  shocks  were  extending  far  enough  so  that  the  stiff  rubber  bump 
stops    would    be    activated. 

B.     TEST  COURSES 

Several    of    the    standard    APG    test     courses     were     used     for     the 
instrumented    test.  These    courses    are    described    below. 

1.        Belgian    Block    Course 

This  course  is  designed  using  "Belgian  Blocks"  to  produce  a 
road  that  has  high  frequency  random  bumps.  The  magnitude  of  each 
bump  is  small.  This  course  is  shown  in  Figure  6-4.  During  the  test 
this   course   was   traversed   at   speeds   ranging   from    10   to   25   MPH. 
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Figure  6-4     APG   Belgian   Block  Course 

2.        Six    Inch    Washboard    Course 

This  course   has   six   inch    bumps    that    are    evenly    spaced    six 
feet  or  72   inches  apart  as   shown  in  Figure  6-5.       This    course    produces 
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significant  heave  and  pitch  effects  on  wheeled  vehicles  even  at 
relatively  slow  speeds.  During  the  test  this  course  was  traversed  at 
speeds   ranging    from    three   to   six    MPH. 
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Figure   6-5      APG   Six    Inch   Washboard   Course 

3.  Cross    Country    Course    Number    1 

This  course  is  similar  to  a  relatively  smooth  but  winding 
trail    road.      This   course   was   traversed   at   approximately   25    MPH. 

4.  Cross    Country    Course    Number    3 

This  is  similar  to  a  rough  trail  road  which  has  random  large, 
medium,  and  small  sized  bumps.  This  course  was  traversed  at 
speeds   ranging   from   5   to    15   MPH. 
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C.  TEST  PROCEDURES 

Each  of  the  above  courses  were  traversed  by  the  howitzer  being 
towed  by  the  HUMMV  at  speeds  within  the  ranges  listed.  Data  was 
collected  on  all  the  channels  described  above,  using  a  sampling  rate 
of  2,500  Hertz.  Because  of  memory  limitations  in  the  data  collection 
equipment  only  approximately  30  seconds  worth  of  continuous  data 
could  be  stored.  For  each  course  several  data  runs  were  completed 
using  various  speeds.  For  instance,  on  the  Six  Inch  Washboard 
Course,  six  data  runs  were  taken  at  speeds  ranging  from  three  to  six 
MPH. 

The  speeds  for  each  data  run  on  the  Six  Inch  Washboard  Course 
and  Cross  Country  Course  Number  3  were  measured  using  the 
HUMMV's  speedometer,  because  these  courses  were  too  rough  to  use 
more  accurate  devices.  For  the  data  runs,  on  these  courses,  the 
driver  tried  to  hold  the  HUMMV  at  the  specified  speed  while  himself 
being     bounced     around.  Because     of     the     extremely     rough      road 

conditions  and  inaccuracies  in  the  speedometer,  the  actual  speeds 
ended  up  being  different  than  the  speeds  specified  in  the  test  plan. 
Consequently,  the  only  accurate  method  of  determining  the  actual 
speed  during  a  data  run  on  these  courses  is  from  the  video  tape 
taken    during    the    run. 

D.  TEST  RESULTS   AND  DISCUSSION 

The   Belgian    Block   and   Cross    Country    Course    Number     1  showed 

no    significant    levels    of    strain    or    acceleration    even    at    speeds  up    2  5 

MPH.       These    tests    simply    proved     that     the     howitzer     could  traverse 
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these  courses  without  causing  damage  to  the  howitzer.  The  Six  Inch 
Washboard  Course  and  Cross  Country  Course  Number  3,  however,  did 
produce  significant  levels  of  strain  and  acceleration  on  the  howitzer 
components. 

The  test  results  for  the  Six  Inch  Washboard  Course  were 
analyzed  in  both  the  time  and  frequency  domains.  Data  runs  at 
speeds  of  3.75,  5.0  and  5.25  MPH  were  analyzed  for  the  Six  Inch 
Washboard  Course  and  compared  to  the  predicted  results  from  the 
computer  model.  The  most  severe  data  run  from  Cross  Country 
Course  Number  3  at  15  MPH  was  analyzed  in  the  time  domain.  The 
procedures   used   and   the   analysis   of  these   results   are   as   follows: 

1.        Time    History    Results 

The  time  history  for  vertical  "in  plane"  A-Frame  force  in 
each  case  was  produced  by  first  summing  the  "in  plane"  axial  forces 
from  both  of  the  A-Frame  legs.  The  same  calculation  procedures  as 
in  Chapter  V  were  used  for  determination  of  the  axial  forces  in  the 
legs  from  the  strain  gage  data  for  strain  gage  pairs  1  and  4  (see 
Figure  5-6).  The  vertical  force  was  determined  by  multiplying  the 
axial  force  times  the  cosine  of  28°.  This  total  vertical  force  in  the  A- 
Frame  legs  was  multiplied  by  the  calibration  factor  of  0.88  that  was 
determined     in     Chapter     V.  The     A-Frame     lateral      forces      were 

calculated  by  determining  the  difference  in  the  horizontal  force 
components  of  the  axial  force  in  each  A-Frame  leg.  The  accelerations 
were    plotted    directly    from    the    measured    data. 

The  force  and  acceleration  time  histories  for  the  washboard 
course    (Figures    6-6    through    6-11)     show     that    high     frequency     spikes 
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occur  on  the  washboard  course.  These  spikes  are  caused  the  braking 
system  being  activated  as  the  howitzer  went  over  the  bumps  for  all 
speeds.  Also,  for  the  5.25  MPH  data  run  the  spade  hit  the  tops  of  the 
bumps.  These  effects  are  not  included  in  the  computer  model  and 
result  in  differences  between  the  model  results  and  the  actual 
results,  particularly,  for  the  5.25  MPH  data  because  the  spade  was 
hitting. 

Time  histories  for  A-Frame  force  for  speeds  3.75,  5.0  &  5.25 
MPH  from  the  actual  test  (Figures  6-6,  6-8,  &  6-10)  are  similar  to  the 
time  histories  predicted  by  the  computer  model  for  the  same  speeds 
(Figures  4-1,  4-3  &  4-5).  The  high  frequency  force  spikes  do  not 
appear  in  the  model's  predicted  results,  however,  the  levels  of  the 
basic  force  response  in  the  actual  results  closely  compare  to  the  force 
response  of  the  model.  For  example,  the  actual  results  for  5.0  MPH 
(Figure  6-8)  has  spikes  that  reach  a  maximum  force  of  5,500  lbs.  But 
the  basic  force  response  on  average  reaches  approximately  2,400  lbs 
in  tension  and  4,300  lbs  in  compression.  The  corresponding  model 
prediction  (Figure  4-3)  shows  forces  of  2,800  lbs  in  tension  and 
4,800  lbs  in  compression.  A  similar  comparison  can  be  made  at  the 
3.75  and  5.25  MPH  speeds.  Although  the  actual  results  are  not 
exactly  the  same  as  the  model  predictions,  they  are  close  enough  to 
conclude  that  the  model  is  reasonably  representative  of  the  actual 
results. 

The  A-Frame  vertical  force  results  for  the  most  severe 
portion  of  the  Cross  Country  Course  Number  3  are  shown  in  Figure  6- 
12.      These  results   show  that    on    a    rough    trail    road    the    A-Frame    force 
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is     not     as     great     as     the     force     that     was     caused    by     the     Six     Inch 
Washboard   Course   at   a   speed  of  5.25    MPH. 
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Vertical  Force  from  APG  Test  (3.75  MPH) 


Figure   6-6      A-Frame   Force  at  3.75   MPH   on   the  Washboard  Course 
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Vertical  Trail  Accel  from  APG  Test  (3.75  MPH) 
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Figure     6-7  Trail     Vertical      Acceleration      at     3.75      MPH      on      the 

Washboard     Course 
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A-Frame  Vertical  Force  from  APG  Test  (5.0  MPH) 
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Figure   6-8      A-Frame  Force  at  5.0  MPH  on  the   Washboard   Course 
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Vertical  Trail  Accel  from  APG  Test  (5.0  MPH) 
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Figure  6-9     Trail  Vertical    Acceleration    at    5.0    MPH    on    the    Washboard 
Course 
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A-Frame  Vertical  Force  from  APG  Test  (5.25  MPH) 
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Figure   6-10      A-Frame   Force   at   5.25   MPH   on   the   Washboard   Course 
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Figure      6-11      Trail     Vertical      Acceleration      at      5.25      MPH      on      the 
Washboard     Course 
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Vertical  Force  for  Cross  Coutry  #3(15  MPH) 


Figure    6-12        A-Frame    Force    at    15    MPH     on     Cross    Country     Course 
Number    3 
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2.         Frequency     Domain 

The  frequency  domain  analysis  for  the  actual  A-Frame 
force  and  trail  acceleration  results  are  shown  in  Figures  6-13  through 
6-18.  The  same  procedures  using  the  "FFT"  function  in  Matlab  as 
were  used  in  Chapter  IV  for  the  model  results  were  used  to  perform 
the  frequency  analysis  on  the  actual  results.  The  time  increment 
used  was  1/1250  seconds  and  was  determined  by  using  every  other 
data  point  of  the  original  data  that  was  taken  at  2500  Hz.  8192 
points  were  used  to  achieve  an  adequate  frequency  discretization. 
The  record  length  was  8192/1250  and  frequency  increment  was 
1250/8192.  The     FFT     analysis     of     both     A-Frame     force     and     trail 

vertical  acceleration  for  each  data  run  clearly  shows  two  frequency 
spikes.  The  lowest  frequency  and  largest  spike  corresponds  to  the 
forcing  frequency  or  road  input  to  the  howitzer  tires.  The  second 
largest    spike    corresponds    to    twice    the    forcing    frequency. 

The  frequency  analysis  for  actual  A-Frame  force  for  speeds 
of  3.75,  5.0,  &  5.25  MPH  (Figures  6-13,  6-15,  &  6-17)  compares 
reasonably  well  to  the  corresponding  frequency  analysis  from  the 
model  results  (Figures  4-2,  4-4,  &  4-6).  For  example,  the  FFT  results 
from  the  actual  data  at  5.25  MPH  (Figure  6-17)  shows  a  large  force 
component  1180  lbs  at  the  road  input  frequency  (1.28  Hz)  and 
second  large  force  component  of  575  lbs  at  double  the  road  input 
(2.56  Hz).  Where  as  the  FFT  of  the  model  for  the  corresponding 
predicted  A-frame  force  (Figure  4-6)  shows  a  large  force  component 
of   1375   lbs   at   the   road   input   frequency   (1.28   Hz)   and   a   second   large 
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force  component  of  500  lbs  at  twice  the  input  frequency  (2.56  Hz). 
For  this  case,  the  total  amount  of  force  from  the  1.28  and  2.56  Hz 
spikes  for  the  actual  results  and  the  model  results  are  approximately 
the  same  (1765  lbs  versus  1875  lbs).  Although  FFT  results  for  the 
actual  data  and  the  model  predictions  are  not  exactly  the  same,  they 
are    close    enough    further    to    validate    the    model's    predicted    results. 
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Figure  6-13      FFT  of    A-Frame    Force    at    3.75    MPH    on    the    Washboard 
Course 
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Vertical  Trail  Accel  FFT  from  APG  Test  (3.75  MPH) 
0.3 


0.25 


0.2 


(A 

o 

V     0.15 

8 
< 


0.1 


0.05 


0 


i i ^*  i jl. i 


0 


Freq  (Hz) 


10 


Figure    6-14       FFT    of   Trail    Vertical    Acceleration    at    3.75    MPH    on    the 
Washboard     Course 
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Vertical  Force  FFT  from  APG  Test  (5.0  MPH) 
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Figure  6-15      FFT  of    A-Frame    Force    at    5.00    MPH    on    the    Washboard 
Course 
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Vertical  Trail  Accel  FFT  from  APG  Test  (5.0  MPH) 
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Figure    6-16      FFT  of  Trail    Vertical     Acceleration    at    5.0    MPH    on    the 
Washboard     Course 
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Vertical  Force  FFT  from  APG  Test(5.25  MPH) 
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Figure   6-17      FFT  of    A-Frame    Force    at    5.25    MPH    on    the    Washboard 
Course 
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Figure    6-18    FFT    of   Trail    Vertical    Acceleration    at    5.25    MPH    on     the 
Washboard     Course 
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3.        Other    Results 

The  strain  reading  from  the  strain  rosettes  that  were 
attached  to  the  bottom  the  howitzer's  trail  showed  low  strain 
readings  for  all  courses.  These  low  strain  readings  clearly  showed 
that  transporting  the  howitzer  in  the  firing  position  with  the  A- 
Frame  attached  produces  insignificant  levels  of  stress  in  the  howitzer 
trail. 

Lateral  forces  in  the  A-Frame  were  also  determined  from 
the  actual  test  data.  Figures  6-19  through  6-21  show  the  lateral 
force  results  for  3.75,  5.0  &  5.25  MPH,  respectively.  The  maximum 
lateral  forces  were  less  than  ten  percent  of  the  maximum  vertical 
forces  for  all  runs  on  the  washboard  course.  For  the  rough  terrain 
Cross  Country  Course  Number  Three  the  maximum  lateral  force 
(Figure  6-22)  was  twenty  percent  of  the  maximum  vertical  force. 
These  relatively  low  values  of  lateral  force  help  to  show  the  validity 
of  the  modeling  technique  used,  since  the  computer  model  only 
predicts  the  vertical  forces  in  the  A-Frame  and  assumes  that  the 
lateral    forces    would    be    comparatively    low. 
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A-Frame  Lateral  Force  from  APG  Test  (3.75  MPH) 


Figure    6-19      Lateral    Force    from   the    APG    Washboard    Course    at    3.75 
MPH 
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A-Frame  Lateral  Force  from  APG  Test(5.0  MPH) 
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Figure    6-20       Lateral    Force    from    the    APG    Washboard    Course    at    5.0 
MPH 
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A-Frame  Lateral  Force  from  APG  Test  (5.25  MPH) 
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Figure    6-21       Lateral    Force    from   the    APG    Washboard    Course    at    5.25 
MPH 
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A-Frame  Lateral  Force  for  Cross  Coutry  #3(15  MPH) 


Figure    6-22       Lateral    Force    from    Cross    Country    Course    Number    Three 
at    15   MPH 
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VII.        CONCLUSIONS/RECOMMENDATIONS 

A.      CONCLUSIONS 

The  computer  model  developed,  produced  output  results  that  are 
reasonably  representative  of  the  actual  results  for  the  washboard 
type  road  input  used  at  Aberdeen  Proving  Grounds.  Therefore,  the 
model  is  validated  and  can  be  used  to  predict  "worst  case"  loading  on 
the  howitzer's  critical  parts.  These  "worst  case"  predictions  for  a 
trench  type  pothole  show  the  Ml  19  Howitzer  can  be  safely  towed  in 
the    firing    position    with    the    A-Frame    attached. 

As  an  additional  result  of  this  study,  we  believe  that  the  T-Bar  is 
the  most  critical  component  on  the  howitzer  in  regards  to  mechanical 
failure  when  towing  in  the  firing  position.  The  predicted  factor  of 
safety  for  the  T-bar  based  on  the  "worst  case"  pothole  terrain  is  1.81. 
Due  to  the  high  cost  of  repair,  making  the  T-Bar  out  of  thicker  steel 
would    provide    an    increased    factor   of   safely. 

The  predicted  factor  of  safety  for  the  A-Frame,  from  the  same 
"worst  case"  terrain,  is  1.92.  This  indicates  that  the  A-Frame  is 
strong  enough  to  withstand  the  worst  case  expected  loads.  However, 
the  A-Frame  is  made  out  of  aluminum  and  aluminum  does  not  have 
a  fatigue  limit.  Thus,  eventually  the  A-Frame  is  expected  to  fatigue. 
The  A-Frame's  fatigue  life  could  be  greatly  extended  by  eliminating 
the  holes   in   its   legs. 
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The  strain  reading  taken  from  the  trail  rosettes  during  the 
instrumented  test  at  APG  show  that  the  stress  in  the  trail  is 
insignificant    during    towing    in    the    A-Frame    position. 

B.      RECOMMENDATIONS 

The  US  Army  should  remove  the  towing  restrictions  and  allow 
the  howitzer  to  be  towed  in  the  firing  configuration  with  the  A- 
Frame    attached    during    tactical    operations. 

The  T-Bar  should  be   made  out    of    thicker    steel    to    provide    added 
assurance   that   it   will    not   fail. 

The  A-Frame  should  be  re-configured  for  manufacture  without 
the  hole  in  its  legs. 
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APPENDIX       A 
SYSTEM    PARAMETER    VALUES 

The  following  table  defines  the  HUMMV/Howitzer  system 
parameters  using  an  inches  (in),  pounds  (lbs),  and  seconds  (sec) 
system  of  units.  Gravity  (G)  is  approximated  to  be  equal  to  386.4 
in/ sec2  .  The  nomenclature  for  the  parameters  that  are  identified  in 
this  Appendix  is  the  same  nomenclature  that  is  used  throughout  this 
paper.  Refer  to  Figures  2-1  and  2-2  to  identify  the  parameters 
defined  below.  The  values  listed  in  the  table  incorporate  the  effects 
of  both  sides  for  those  components  that  are  identical  on  the  right  and 
left  side  of  the  HUMMV  or  Howitzer.  For  example,  the  spring 
stiffness  for  each  of  the  front  HUMMV  tires,  Kx,  is  1323.0  lbs/in.  This 
value  is  multiplied  by  two  to  account  for  the  combined  effect  of  both 
front  HUMMV  front  tires.  The  value  listed  in  the  following  table  for 
Kj  is    2646    lbs/in. 
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TABLE  A-l     COMPUTER  MODEL  PARAMETER  VALUES 


Symbo 

Value/units 

How    Determined/Notes 

Ki 

2646.0  lbs/in 

Reference  1/HUMMV  front  tires  stiffness 

K2 

1908.0  lbs/in 

Reference   1  /HUMMV  front  suspension  stiffness 

*S 

3544.0  lbs/in 

Reference  1/  HUMMV  rear  tires  stiffness 

K4 

5040.0  lbs/in 

Reference  1  /  HUMMV  rear  suspension  stiffness 

Ke 

2000.0  lbs/in 

Estimated /Howitzer  tires  stiffness 

K7 

1056.0  lbs/in 

Reference  2 /Modeled  Howitzer  vertical 
suspension  stiffness 

Ka 

200,000.0  lbs/in 

Calculated/Combine  A-Frame    &  T-Bar  stiffness 

Bt 

Bic=  171.6  Ibs-sec/in 
Bit=308.4  Ibs-sec/in 

Reference  1/Bic  is  for  compression  &  Bit 
for  tension.    Values  are  20%  more  than  Ref  1 

B2 

B2C=  199.2  Ibs-sec/in 
B2T=480.0  Ibs-sec/in 

Reference  1/  B2C  is  for  compression  &  B2T 
for  tension.    Values  are  20%  more  than  Ref  1 

B3 

B3c=40.0  Ibs-sec/in 
B3T=  132.0  Ibs-sec/in 

Calculated  from  RIA  data/  B3C  is  for 
compression  &  B3T  for  tension. 

B4 

250,000  lbs-in-sec/rad 

Estimated/  Torsional  Damping  in  elevating  clutch 

M, 

1.294  Ibm 

Kererence  1/  includes  mass  or  HUMMVtront 
wheels.shocks  and  susDension  sDrinas 

M2 

1.294  Ibm 

Reference  1  /  Includes  mass  of  HUMMV  rear 
wheels.shocks  and  suspension  springs 

M3 

21.22  Ibm 

Reference  1/  HUMMV  Body 

M4 

1.035  Ibm 

Reference  2/lncludes  howitzer  wheels, 
shocks,  and  30%  of  the  road  arm  masses 

M5 

9.834  Ibm 

Reference  2/  Howitzer  trail  &  saddle  mass 

Me 

5.952  Ibm 

Reference  2/  Howitzer  elevating  mass 

h 

52,680.0  lbs-in-sec2 

Reference  1/  Moment  of  inertia  about  HUMMV  CG 

99 


TABLE  A-l     COMPUTER  MODEL  PARAMETER  VALUES  (CONT) 


Symbo 

Value/units 

How    Determined/Notes 

'2 

10,000  Ibs-in-  sec2 

Estimated/  Moment  of  inertia  of  Trail  &  saddle 
mass  about  CG 

'3 

15,000   lbs-in-sec2 

Estimated/  Moment  of  inertia  of  Elevating  mass 
about  CG 

Li 

83.8  in 

Reference    1  /HUMMV  CG  to  front  springs  &  shocks 

L2 

46.7  in 

Reference    1  /HUMMV  CG  to  rear  springs  &  shocks 

L3 

75.7  in 

Reference   1  /HUMMV  CG  to  hitch  connecting  point 

L7 

130.5  in 

Measured/  HUMMV  front  axle  to  HUMMV  rear  axle 

L8 

294.5  in 

Measured/  HUMMV  front  axle  to  Howitzer  axle 

Mi 

104.4  in 

Measured/  Howitzer  axle  to  CG  of  Trail  &  Saddle  mass 

Liz 

135.0  in 

Measured/Howitzer  axle  to  hitch  connecting  point 

Li  3 

38.0  in 

Measured/Trunion  to  where  A-Frame 
connects  to  the  T-Bar 

hi 

110.3  in 

Measured/Trunion  to  hitch  connecting  point 

LA 

45.0  in 

Measured/  Trunion  to  where  A-Frame 
attaches  to  the  Elevating  mass/Gun  Tube 

Lcg 

32.5  in 

Provide  by  RIA/  Trunion  to  CG  of  Elevating  mass 

e0 

1 2.0  degrees 

Measured/  Static  angle  from  horizontal  to 
Gun  Tube  when  A-Frame  is  Attached 

*A 

51.0  degrees 

Measured/  Angle  from  Gun  Tube  to  A-Frame 
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APPENDIX      B 

MATRICES   FOR   EQUATIONS   OF   MOTION 

This  Appendix  contains  the  mass,  damping,  and  stiffness 
matrices  for  the  equations  of  motion  that  were  developed  in  Chapter 
II.      The  mass  system  matrix,   [M],  is  as  follows: 


rmu  0 

0 

0 

0 

0 

o      1 

'0     m22 

0 

0 

0 

0 

0 

'o     0 

m33 

0 

m35 

m36 

m37 

0      0 

0 

m44 

0 

0 

0 

,0      0 

m53 

0 

m55 

m56 

m57 

0      0 

m63 

0 

m65 

m66 

m67 

0      0 

m73 

0 

m75 

m76 

m77    _ 

where: 

mn  =M, 

m22  = 

M2 

m35  =  M5L3    + 

M6L3 

m37  =  M6LCG  cos(0o) 
m53  =  M5L3  +  M6L3 

m56  =M5L3L11     +    M6L3L16 

m63  =M5LU  +  M6L16 

m66  =  I2  +  M5Ln2  +  M6L162 

m73   =  M6LCG  COS(0O) 

m76  =  M6LCGL16cos(80) 


m33  =  M3  +  M5+M6 

m36=M5Lll     +     M6L16 

m44  =  M4 

m55  =  Ix  +   M5L32  +  M6L32 
m57=M6LCGL3cos(0o) 

m65=M5L3Lll    +    M6L3Ll6 
m67  =  M6LCGL16COS(0o) 
m75  =  M6LCGL3COS(eo) 

m77  =  I3  +  M6LCG2 


Note    due     to     symmetry     of     [M]     the     following     matrix     elements     are 


equal 


m53  =  m35 

m65  =  m56 


m63   =  m36 

m75  -  m57 


m73  =  m37 
m76  =  m67 
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The  damping  matrix  for  the  system  is  dependent  on  whether 
the  shock  absorber  forces  are  tensile  or  compressive  at  a  given  time. 
The  tensile  or  compressive  state  for  each  of  the  shock  absorbers  is 
determined  by  the  relative  velocities  of  the  end  connection  points  for 
the   shock   absorber.      The   damping    matrix,    [C(x)],  is  as  follows: 


|Cn 

0 

C13 

0 

lo 

v^T 

C23 

0 

1    31 

C32 

C33 

c3 

0 

0 

C43 

c4 

|C5i 

C52 

C53 

c5 

0 

0 

C63 

c6 

0 

0 

0 

1 

where: 

Ci,  = 

B, 

C22  = 

B2 

c31  - 

-B, 

C34  = 

"B3 

C36  = 

B3L12 

c45  = 

~~  BjLj 

c52  = 

-B2L2 

C54  : 

=  ~B3L3 

C56  = 

B3L3L1: 

» 

c65  = 

BjLjL^ 

C76 

"B4 

34 


54 


'15 


'25 


'35 


'4  5 


'55 


'65 

0 


0 

0  1 

0 

0 

C36 

0 

C46 

0 

C56 

0 

c66 

c67 

C76 

C77_ 

Cis-'B, 

^23  "2 


C32         r>2 


'3  5 


-  -BjLj  +  B2L2  +   B3L3 


C43  =  "B3 

C46  =  ~B3Lj2 


C15  ~  B^! 

c25  =  "B2L2 

C33  -  Bi  +  B2  +  B3 


C44  -  B3 


'51 


B.L, 


c53  ~  ~B1L1  +  B2L2  +    B3L3 


c55  =  BjV  +  B2L/  +  B,V 


'63 


~  BjLi: 


c66  ~    BjL12    +  B4 

C77  -  B4 


c64         B3L12 


c67  =  "B4 


Note   due  to   symmetry   of  [C(x)]    the    following    elements   are    equal: 


C31  ~  C13 
C52  ~~  C25 
C64  =  C46 


C32=C 


23 


C53        C35 

C65  ~~  C56 


C51        C15 


C54        C45 

C76  ~  C67 


C43         C34 


C63        C36 
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The   system   stiffness   matrix,    [K],   is  as   follows: 


0 
0 
0 

Lo 


22 


^32 

0 

0 
0 


k13 
k23 

k33 

^43 

k53 

^63 

0 


0 
0 

k34 
k44 

k54 

K< 
0 


*N5 
^25 
^35 

k45 

k65 
0 


36 

k 


46 


^56 
k66 

k. 


0 

0 
0 

0 
0 

k 


67 


76 


k77  J 


where: 

ku  =  Kt  +  K2 

k22  =  ^3     +      ^4 

kj,  =  -K2 
kj4  =  -K7 
kj6  =  K7L12 
k45  =  -K7L3 
kjj  =  ~K4L2 

k54   =    _K-7L3 
k56  ~~  *^7^3^12 
*65   ~~  ^7^3^12 

k67  =  "  KALAL13cos(0o) 
k77=    KALAL13cos(60) 


kn  =  ~K2 

^23  ^-4 

k32   =  ~^4 


ki5  =  ^2^1 

k15  =  ~K4L2 

k33  =  K2  +  K4  +K7 


k35  -     K2L,    +   K4L2  +   K7L3 

k43  =  ~K7  k44  =  K6  +  K7 

k46  =  -K7L12  k51  -  K2Lj 

53   =  _      2     1      "^  4     2  7     3 

k5S  -  K2L!     +   K4L2    +  K7L5 

k63   ~~  ^7^12  k64  ^7^12 

k66  =  K7L122  +  KALAL13cos(0o) 
k76=  -  KALAL13cos(0o) 


Note   due   to    symmetry   of   [K]    the   following    matrix   elements   are   equal: 


*31  =   ^3 

ks2  ~  k2s 
k64  =  k46 


k32         ^23 

k     =  k 
k     =  k 

*"65         K56 


^si     k15 

k54  =  k45 
k     =  k 

K76        K67 


k43     kj4 
k63  =  k36 
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APPENDIX      C 
NUMERICAL    INTEGRATION    COMPUTER    CODE 

This    Appendix    contains    the    computer    code    used    to    integrate 
the   equations   of   motion    using    the   ODE   45    integration    routine    in 
Matlab.. 

1.      THE   ODE   45   INTEGRATION   FUNCTION 

function     xdot=mateqn(t,x) 

%  Define    gravity 
G=386.4; 

%  Define    HUMMV    Parameters: 

Ll=83.8;  L2=46.7;  L3=75.7;      L7=130.5;      L8=294.5; 

Kl=2*1323.0;        K2=2*954.0;        K3=2*1777.0;      K4=2*2520.0; 

K5=1000000.0; 

B1C=2*1.2*71.5;        B1T=2*1.2*128.25;      B2C=2*1.2*83.0; 
B2T=2*1. 2*200.0; 

Wl=2*250.0;         W2=2*250.0;  W3=8200.0; 

M1=W1/G;  M2=W2/G;         M3=W3/G; 

11=52680;        RW=5400;  FW=3800; 

%  Define    Gun    Parameters: 

L4=123.4;  L5=11.6;  LA=45.0;  LS=42.6;         LCG=32.5; 

Lll=104.4;         L12=135.0;     L13=38.0;  L14=17.0;     L15=9.70; 

L16=110.3; 

K6=2*  1000.0;  K7=2*528;  KA=200000;      KS=300; 

B3C=2*20.0;  B3T=2*66.0;         B4=.5*250000.0; 

W4=2*200;  W5=1500;  W6=2300;  M4=W4/G;        M5=W5/G; 

M6=W6/G; 

12=10000;        13=15000; 

Q0=12.0*pi/180;  PHIA=51.0*pi/180;        PHIS=14.0*pi/180; 

AW=3800;        TB=88000;        %  Define  Mass  matrice 

%  Define    the    mass    matrix 
M=[zeros(7,7)]; 
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M(1,1)=M1;       M(2,2)=M2;    M(3,3)=M3+M5+M6;    M(3,5)=M5*L3+M6*L3; 

M(3,6)=M5*L1 1+M6*L16;    M(3,7)=M6*LCG*cos(Q0); 

M(4,4)=M4;       M(5,3)=M5*L3+M6*L3;    M(5,5)=I1+M5*L3A2+M6A2; 

M(5,6)=M5*L3*L11+M6*L3*L16;       M(5,7)=M6*L3*LCG*cos(Q0); 

M(6,3)=M5*L11+M6*L16;        M(6,5)=M5*L3*L11+M6*L3*L16; 

M(6,6)=I2+M5*L11A2+M6*L16A2+M6*L14A2; 

M(6,7)=M6*L16*LCG*cos(Q0); 

M(7,3)=M6*LCG*cos(Q0);      M(7,5)=M6*L3*LCG*cos(Q0); 

M(7,6)=M6*L16*LCG*cos(Q0); 

M(7,7)=I3+M6*LCGA2; 


%  Define    Stiffness    Matrix 

K=[zeros(7,7)]; 

K(1,1)=K1+K2;     K(l,3)=-K2;     K(1,5)=K2*L1; 

K(2,2)=K3+K4;     K(2,3)=-K4;     K(2,5)=-K4*L2; 

K(3,l)=-K2;     K(3,2)=-K4;     K(3,3)=K2+K4+K7; 

K2*L1  +K4*L2+K7*L3;      K(3,6)=K7*L12; 

K(4,3)=-K7;     K(4,4)=K6+K7;     K(4,5)=-K7*L3; 


K(3,4)=-K7;     K(3,5)=- 
K(4,6)=-K7*L12; 


K(5,1)=K2*L1;      K(5,2)=-K4*L2;   K(5,3)=-K2*L1+K4*L2+K7*L3;      K(5,4)=- 

K7*L3; 

K(5,5)=K2*L1A2+K4*L2A2+K7*L3A2;      K(5,6)=K7*L3*L12; 

K(6,3)=K7*L12;      K(6,4)=-K7*L12;      K(6,5)=K7*L3*L12; 

K(6,6)=K7*L12A2+KA*LA*L13*cos(Q0); 

K(6,7)=-KA*L13*LA*cos(Q0); 

K(7,6)=-KA*L13*LA*cos(Q0); 

K(7,7)=+KA*L13*LA*cos(Q0); 

%  Define    and    velocities: 

VHF=x(10)   -  (Ll*x(12));  VHB=x(10)   +  (L2*x(12)); 

VGB=x(10)    +   (L3*x(12))    +   (L12*x(12)); 

%  Define    Road    Inputs: 

VMPH=5.0;      VEL= VMPH*  1 7.6;      PD=72/VEL; 

Ul=3*(sin(2*pi*t/PD)); 

if  t  <  L7/VEL,     U2=0;else,       T2=t-(L7/VEL); 
U2=3*(sin(2*pi*T2/PD));    end 

if  t  <  L8/VEL,     U3=0;else,       T3=t-(L8/VEL); 
U3=3*(sin(2*pi*T3/PD));    end 
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%    Determine    the    shock    absorber    damping    constants    based    on 
compression     or    tension 

%   HUMMV  Front  shock 
if  (VHF  -  x(8))">  0.0,     B1=B1T;  else,  B1=B1C;  end 

%   HUMMV  Rear  shock 
if  (VHB  -  x(9))  >  0.0,  B2=B2T;  else,     B2=B2C;  end 

%  GUN    shock 
if  (VGB  -  x(ll))  >  0.0,     B3=B3T;  else,     B3=B3C;  end 

%  Define    Damping    Matrix 

B=[zeros(7,7)]; 

B(1,1)=B1;        B(1,3)=-B1;   B(1,5)=B1*L1; 

B(2,2)=B2;       B(2,3)=-B2;       B(2,5)=-B2*L2; 

B(3,1)=-B1;     B(3,2)=-B2;     B(3,3)=B1+B2+B3;     B(3,4)=-B3;     B(3,5)=- 

B1*L1+B2*L2+B3*L3;      B(3,6)=B3*L12; 

B(4,3)=-B3;      B(4,4)=B3;      B(4,5)=-B3*L3;      B(4,6)=-B3*L12; 

B(5,1)=B1*L1;      B(5,2)=-B2*L2;      B(5,3)=-B1*L1+B2*L2+B3*L3; 

B(5,4)=-B3*L3;      B(5,5)=B1*L1A2+B2*L2A2+B3*L3A2; 

B(5,6)=B3*L3*L12; 

B(6,3)=B3*L12;      B(6,4)=-B3*L12;      B(6,5)=B3*L3*L12; 

B(6,6)=B3*L12A2+B4; 

B(6,7)=-B4;       B(7,6)=-B4;       B(7,7)=B4; 

%  Define    Forcing    Vector 
F=[zeros(7,l)]; 

F(1,1)=K1*U1;      F(2,1)=K3*U2;   F(4,1)=K6*U3; 

%  CHECK  IF  FRONT  HUMMV  WHEEL  LOSES  CONTACT 
if  K1*(U1  -  x(l))  <  -FW,       F(2,1)=-FW;       K(1,1)=K2;  end 

%  CHECK  IF  REAR  HUMMV  WHEEL  LOSES  CONTACT 

if  K3*(U2  -  x(2))  <  -RW,       F(2,1)=-RW;       K(2,2)=K4;  end 

%  CHECK  IF  HOWITZER  WHEEL  LOSES  CONTACT 

if  K6*(U3  -  x(4))  <  -AW,       F(4,1)=-AW;       K(4,4)=K7;  end 

%  transform    into    system    of    1st    order    ode's 

id=[l    zeros(l,6);0    1    zeros(l,5);0   0    1    zeros(l,4);zeros(l,3)    1    zeros(l,3)]; 

id=[id;zeros(l,4)    1    0   0;zeros(l,5)    1    0;zeros(l,6)    1]; 
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a=[zeros(7,7)    id;    -inv(M)*K    -inv(M)*B]; 
fl=inv(M)*    F; 
f=[[zeros(7,l)];    fl]; 

^integrate 
xdot=a*x  +  f; 


2.     THE  DRIVER  CODE  FOR  ODE  45 

%    Set  the   Initial   Condition 
x01=[zeros(14,l)]; 


%    Integrate    in    one    second    increments 
t01=0.00;    tf  1=0.999;     %x01=xn(length(xn),:); 
[tl,xl]=ode45('mateqn\t01,tfl,x01,le-3); 


102=1.0001;  tf2=1.999 
[t2,x2]=ode45('mateqn 
t03=2.0001;  tf3=2.999 
[t3,x3]=ode45('mateqn 
t04=3.0001;  tf4=3.999 
[t4,x4]=ode45('mateqn 
t05=4.0001;  tf5=4.999 
[t5,x5]=ode45('mateqn 


x02=xl(length(xl),:) 
',t02,tf2,x02,le-3); 

x03=x2(length(x2),:) 
\t03,tf3,x03,le-3); 

x04=x3(length(x3),:) 
■,t04,tf4,x04,le-3); 

x05=x4(length(x4),:) 
\t05,tf5,x05,le-3); 

x06=x5(length(x5),:) 


t06=5.0001;     tf6=5.999 
[t6,x6]=ode45('mateqn',t06,tf6,x06,le-3); 
%t07=6.0;    tf7=6.999;    x07=x6(length(x6),:); 
%[t7,x7]=ode45('mateqn\t07,tf7,x07,le-3); 
%t08=7.0;    tf8=7.999;    x08=x7(length(x7),:); 
%[t8,x8]=ode45(,mateqn,,t08,tf8,x08,le-3); 


%    Combine    the    solution    state    matrix 
xn=[xl;x2;x3;x4;x5;x6];%x7;x8]; 
tn=[tl;t2;t3;t4;t5;t6];%t7;t8]; 
save   ghl    xn   tn 
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APPENDIX    D 
POTHOLE  COURSE  COMPUTER  CODE 

This    Appendix    contains    the    computer    code    that    was    developed 
to   model    the   terrain   inputs   for   a   road   with   a   trench   type   pothole.    This 
program    calculates   and   plots   a   trace   of  a   pothole   profile   and   the   axle 
path    followed    by    circular    wheel    rolling    through    this    type    pothole.    This 
program    works    for   any   pothole   that   can    be    model    with    a   flat   bottom, 
and   sides   that   can    be   approximated   as    straight   lines.      The   portion   of 
the    code    shown    below    only    works    for   potholes    that    have   an    entrance 
height   that   is   equal    to   the   exit   height.      Additional   code   to   determine 
the   terrain   inputs   for  a  any   type  of  pothole   is  on   file  at   the  Naval 
Postgraduate    School.       Computer   code    was    also    developed    to    model 
bump   inputs,   such   a   log   in   a  road  and   maintained  on   file   at   the  Naval 
Postgraduate     School. 

COMPUTER  CODE: 

function      [x,y,ya,u]=pothole(xent, depth, thetal,xbot,theta2,farht,xtrail); 

%      The   call   statement   for  this   program/function   is   as    %  follows: 
%  [x,y,ya,u]=pothole(xent, depth, thetal,xbot,theta2,farht,xtrail) 

%  Input   parameters   used   in   the   call   are   defined   as   follows: 

%  depth      -   the  vertical   distance  from  the  bottom  to  the  level 

%  entrance     surface. 

%  thetal    -  the  angle  of  the  decline  on  the  entrance  side,   measured 

%  from     an  extension  of  the  level  entrance  surface  to  the  decline. 

%  xbot       -  the  horizontal  distance  across  the  bottom. 

%  theta2   -   the  angle  decline  on  the  exit  side,   measured  like  thetal. 

%  farht     -  the  vertical   distance  from  the  bottom  to  the  level  exit 

%  surface. 

%  xent        -  the  horizontal  length  of  the  level  entrance  region. 

%  xtrail  -  the  horizontal  length  of  the  level  exit  region. 

%  Output  parameters 

%  x     -  a  column  vector  containing  a  discrete  amount  of  evenly   spaced  'x' 

%  values  that  were  used  to  calculate  the  profile  and  axle  path. 

%  y     -  a  column  vector  that  contains  the  'y'  values  for  the  pothole 
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%  profile.     The  bottom  of  the  pothole  is  defined  as  "zero"  y  value. 

%  ya  -   a  column   vector  containing   the  axle  path   'y'   values. 

%  u     -  a  column  vector  containing  the  height  of  the  bottom  of  the 
%  wheel.  V  is  the  input  for  the  simulation  model 

%  Other    parameters    assigned    in    the    program: 

%        deltax  -  the  V  intraval   spacing. 

%        xent        -  the  horizontal  length  of  the  level  entrance  region. 

%        xtrail  -  the  horizontal   length  of  the  level  exit  region. 

%        rad  -  the  radius  of  the  wheel  that  is  rolling  through  the  pothole. 

%  Varibles   that   are   calculated   based   on   the   above   assigned 

%  parameters: 

% 

%  rhtetal    &  rtheta2   -   the  thetal    &   theta2   in  radians. 

%  xnear  &  xfar  -  the  x  distance  differance  from  the  top  to  the 

%  bottom  for  the  near  &  far  sides. 

%  xcritl   &  xcrit2  -  the  x  value  where  the  wheel  can  no  longer 

%  follow  a  circular  path  as  it  goes  around  the  near  &  far 

%  corners(near  into  hole,  far  out  of  hole). 

%  xcrl    -  the  x  value  where  the  wheel  hits  the  far  side  or  corner 

%  xcr2  -  the  x  value  from  where  the  wheel  hits  the  the  exit  region 

%  from  the  near  side  or  near  corner. 

%  xchl   -  the  x  value  at  the  bottom  corner  on  the  entrance  side. 

%  xch2  -  the  x  value  at  the  bottom  corner  on  the  exit  side. 

%  xch3  -  the  x  value  at  the  top  corner  on  the  exit  side. 

%  xend  -  the  ending  x  value  of  the  pothole  profile. 

%  ycornl   &  ycorn2  -  the  y  values  when  the  axle  is  following  a 

%  circular  path  around  the  near  &  far  corners  respectively 

%********assignment   of  pothole   &   wheel   parameters****************** 
deltax=.01;  rad=17.0; 

rtheta  1  =theta  1  *pi/ 1 80;,  rtheta2=theta2*pi/l  80; 

xnear=depth/tan(rtheta  1 ) ;,  xfar=farht/tan(rtheta2) ; 

xchl=xnear+xent;,  xch2=xchl+xbot; 

xch3=xch2+xfar;,  xend=xch3+xtrail; 
xcritl  =xent+rad*sin(rthetal);,  xcrit2=xch3-rad*sin(rtheta2); 

xcrit3=xent+rad*sin(rtheta2);,  xcrit4=xch3-rad*sin(rtheta2); 

%    converting/rounding    x-switch    values    to    integer    increments    of 
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%    deltax;    keeps   the   output   x,   y,&ya   vectors   the  same   size 

num=xent/deltax;,   num=round(num);,  xent=num*deltax; 

num  l=xch  l/deltax;,num  l=round(num  1);,  xch  l=num  1  *deltax; 

num2=xch2/deltax;,num2=round(num2);,  xch2=num2*deltax; 

num3=xch3/deltax;,num3=round(num3);,  xch3=num3*deltax; 

num4=xend/deltax;,num4=round(num4);,  xend=num4*deltax; 

xdl=xcritl-xent;,nl=xdl/deltax;,nl=round(nl);,  xcritl=xent     + 
nl*deltax; 

xd2=xch3-xcrit2;,n2=xd2/deltax;,n2=round(n2);,  xcrit2=xch3     - 
n2*deltax; 

or  ****************  POTHOLE  PROFILE  *************************** 

%   ******  creating  x  vectors  for  pothole  regions  ***************** 

xl=0:deltax:xent;,  sizl=size(xl); 

x2=xent+deltax:deltax:xchl;, 

x3=xchl+deltax:deltax:xch2;,  siz3=size(x3); 

x4=xch2  + deltax:  deltax  :xch3;, 

x5=xch3+deltax:deltax:xend;,  siz5=size(x5); 

x=[xl   x2  x3  x4  x5]; 

%  ******  calculating  y  values  for  the  pothole  profile  **************** 

yl=depth*ones(l:sizl(2)); 

y2=     -tan(rthetal)*x2     +tan(rthetal)*xent     +     depth; 

y3=zeros(l:siz3(2)); 

y4=    tan(rtheta2)*x4    -    tan(rtheta2)*xch2    ; 

y5=     farht*ones(l:siz5(2)); 

y=[yl    y2   y3   y4   y5]; 

Of  ****************    AYT  F  PATH  ************************************** 

%******calculating    the   axle   trace   in   the   entrance   regi0n******** 
xal=0:deltax:xent;,sizal=size(xal); 
yal=(depth+rad)*ones(l:sizal(2)); 
if  xcrit2  <=  xcritl, 

%***checking  if  the  axle  motion  is  only  circular***** 
if  xcrit2  >  xent, 

%  finding  the  x  intercept  of  the  near  &  far  side  circular  motion 
for  xint=xcrit2: deltax :xch3, 

yl=depth  +  sqrt((rad.A2)  -  ((xint-xent).A2)); 
y2=farht  +  sqrt((rad.A2)  -  ((xint-xch3).A2)); 
if  yl<=y2;  break,end 
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else 


end 
else 

for  xint=xent+deltax:deltax:xch3, 

yl=depth  +  sqrt((rad.A2)  -  ((xint-xent).A2)); 
y2=farht  +  sqrt((rad.A2)  -  ((xint-xch3).A2)); 
if  yl<=y2;  break,end 
end 
end 

%  calculating  the  near  and  far  side  y  values  for  the  axle  path 
xa2=xent+deltax:deltax:xint; 
ya2=depth  +  sqrt((rad.A2)  -((xa2-xent).A2)); 
xa3=xint+deltax:deltax:xch3; 
ya3=farht  +  sqrt((rad.A2)  -((xa3-xch3).A2)); 
yac=[ya2    ya3]    ; 

xa4=xch3+deltax:deltax:xend;,       siza4=size(xa4); 
ya4=(farht+rad)*ones(l:siza4(2)); 
ya=[yal  yac  ya4]; 
%  checking  &  calculating  axle  path  when  its  not  just  circular 

%****intializes  a  counter  &  the  vector  yac**** 
n=0;,         nnum=(xch3-xent)/deltax;,  yac=zeros(l,nnum); 

%*****checks  &  determination  of  axle  path**** 
for  xa=xent+deltax:deltax:xch3, 
if  xa  <=  xcritl, 

ycornl=depth  +  sqrt((rad.A2)  -  ((xa-xent).A2)); 

yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2+rad/cos(rtheta2); 
if  ycornl  >  rad    &    ycornl  >  yfar,  yacl=ycornl; 
elseif  yfar  >  rad,  yacl=yfar; 
else  yacl=rad;    end 
elseif  xa  >=  xcrit2, 

ycora2=farht  +  sqrt((rad.A2)  -  ((xa-xch3).A2)); 
ynear=-tan(rthetal)*xa+tan(rthetal)*xent... 
+depth +rad/cos(rtheta  1 ) ; 
if  ycorn2  >  rad    &    ycorn2  >  ynear,  yacl=ycorn2; 
elseif  ynear  >  rad,    yacl=ynear; 
else     yacl=rad;    end 


else 


ynear=-tan(rthetal)*xa+tan(rthetal)*xent... 
♦depth +rad/cos(rtheta  1 ) ; 
yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2+rad/cos(rtheta2); 
if  ynear  >  yfar    &    ynear  >  rad,    yacl=ynear; 
elseif  yfar  >  rad,    yacl=yfar; 
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else  yacl=rad;     end 
end 
n=n+l;,  yac(l,n)=yacl; 

end 

xa4=xch3+deltax:deltax:xend;,       siza4=size(xa4); 
ya4=(farht+rad)*ones(l:siza4(2)); 
ya=[yal  yac  ya4]; 
end 
%**********  PLOT  OF  POTHOLE  PROFILE  &  AXLE  PATH  ************* 

ax=[0,xend,-50,yplot]; 

y=y-depth;      ya=ya-depth;u=ya-17; 

x=x';y=y,;ya=ya';u=u'; 

matrix=fx   y   ya   u]; 

save    road2dat.m     matrix 

axis(ax); 

plot(x,y);   hold  on    ;   plot(x,ya)    ,pause,hold   off 
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APPENDIX      E 
NATURAL   FREQUENCIES    AND   MODE   SHAPES 

Table  E-l  below  contains  the  natural  frequencies  and  mode 
shapes  of  the  HUMMV/Howitzer  system.  They  were  calculated  using 
the  procedures  decsribed  in  Chapter  IV.  In  Table  E-l  the  system's 
natural  frequencies  are  listed  in  the  top  row.  The  positional 
coordinate  values  of  the  mode  shape  which  corresponds  to  a 
particular  natural  frequency  are  listed  in  the  column  below  the 
natural  frequency.  The  natural  frequencies  shown  in  Table  E-l  have 
units  of  hertz.  The  positional  coordinates  have  units  of  inches  for  Zi, 
2s,Z3,  and  Z4,   and   radians   for   the   rotational   cordinates   0i,  02,  and  03. 


TABLE  E-l      NATURAL  FREQUENCIES  AND  MODE  SHAPES 

Natural     Frequencies 


Coordi- 
nates 

1,286 

1.870 

3.328 

12.26 

13-61 

18.51 

23.80 

Zi 

-.0001 

.0783 

-.2392 

-.0212 

1.217 

.0130 

.0342 

z2 

.0130 

.1293 

.1425 

.0076 

.0364 

-1.222 

-.1221 

z3 

.0140 

.2057 

-.0408 

.0018 

-.0221 

.0138 

.0277 

z4 

-.1064 

.0168 

.0202 

-1.385 

-0.221 

-.0097 

.0341 

6i 

.0002 

.0003 

.0059 

.0001 

.0011 

.0005 

.0025 

e2 

-.0025 

-.0013 

-.0026 

.0001 

-.0003 

-.0001 

-.0036 

e3 

-.0025 

-.0013 

-.0027 

.0003 

-.0006 

-.0011 

.0063 
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APPENDIX       F 
A-FRAME    CALIBRATION    DATA 

The    following    tables    contain    the    load    versus    strain    data    that 
was    recorded    during    the    A-Frame    calibration    procedures    that    are 
described   in   Chapter   V.      The   locations   of  the   strain   gages   is   shown   on 
Figures   5-1    and   5-6.      Tables   F-l    and   F-2   contain   the   data   from   the 
vertical    load   tests   and   Tables   F-3   and   F-4   contain    the   data    from    the 
lateral    load    tests.      The    method    used    to   calculate    stress   values    from 
the      strain   gage   readings   is   described   in    Chapter   V. 
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TABLE  F-l     VERTICAL  TEST  #1     A-FRAME  CALIBRATION  DATA 
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TABLE  F-l  VERTICAL  TEST  #1  A-FRAME  CALIBRATION  DATA  (CONT) 
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TABLE  F-2     VERTICAL  TEST  #2     A-FRAME  CALIBRATION  DATA 
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TABLE  F-2  VERTICAL  TEST  #2  A-FRAME  CALIBRATION  DATA  (CONT) 
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TABLE  F-3     LATERAL  TEST  #1     A-FRAME  CALIBRATION  DATA 
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TABLE  F-3     LATERAL  TEST  #\  A-FRAME  CALIBRATION  DATA  (CONT) 
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APPENDIX      G 

HOWITZER    ORIENTATION    PICTURES    AND    DIAGRAMS 

This    Appendix    contains    pictures    and    diagrams    that    identify    the 
key    components    of   the    Ml  19    Howitzer    that    were    important    for    this 
study. 


Figure   G-l    Towing   in   Firing   Position 
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Figure    G-2   Towing    in    Stowed   Position 
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Figure    G-3       A-Frame/T-Bar    Support 


124 


Figure  G-4      T-Bar  (View   from  top  Down) 
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Figure    G-5       Howitzer   Wheels    and    Suspension 
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